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STUDY OF EFFECTS OF SWEEP ON THE FLUTTER OF CANTILEVER WINGS 1 

By J. G. Barmby, H. J. Cunningham, and I. E. Garrick 


SUMMARY 

An experimental and analytical investigation of the flutter 
of sweptback cantilever wings is reported. The experiments 
employed groups oj wings swept back by rotating and by 
shearing. The angle of sweep ranged from 0 ° to 60° and 
Mach numbers extended to approximately 0.85. A theoretical 
analysis of the air forces on an oscillating swept wing of high 
length-chord ratio is developed, and the approximations in- 
herent in the assumptions are discussed. Comparison with 
experiment indicates that the analysis developed in the present 
report is satisfactory for giving the main effects of sweep , at 
least for nearly uniform cantilever wings of high and moderate 
length-chord ratios. A separation of the effects of finite span 
and compressibility in their relation to sweep has not been made 
experimentally but some combined effects are given. A dis- 
cussion of some of the experimental and theoretical trends is 
given with the aid of several tables and figures. 

INTRODUCTION 

The present report is an outgrowth of the trend toward 
the use of swept wings for high-speed flight and presents the 
results of an analysis and of an accompanying exploratory 
program of research in the Langley 4.5-foot flutter research 
tunnel on swept cantilever wings. The material was 
assembled in a memorandum form with a similar title in 
1948. The chief purposes of the present report are to 
provide a more detailed exposition of the analysis and to 
make the main material more generally available. 

Some previous experimental and analytical work on swept 
wings is mentioned here. A preliminary experimental 
investigation of the effect of sweep on flutter has been made 
(reference 1) with a single, simple rigid wing mounted flexibly 
on a base which could be rotated to various desired sweep 
angles. This investigation was made at low Mach numbers 
for two bending-torsion frequency ratios and at several 
angles of sweepback. Another investigation (data un- 
published) in which the density of the test medium was a 
variable was conducted by D. Benun on the same type of 
rigid, flexibly mounted wing at higher Mach numbers and at 
sweep angles of 0° and 45°. Other unpublished work on 
swept wings exists, but a search of the available information 
indicates a need for further systematic study. 


The experimental work reported herein dealt with models 
mounted as cantilevers at their roots. These cantilever 
models differed from the rigid, flexibly mounted wings, 
which had all bending and torsion flexibility concentrated 
at the root, and thus were subject to different root effects. 
In order to facilitate analysis the cantilever models w'ere 
uniform and untapered. The intent of the experimental 
program was to establish trends and to indicate orders of 
magnitude of the various effects of sweep on flutter rather 
than to isolate precisely^ the separate effects. 

The models were swept back in two basic manners — 
shearing and rotating. For the case in which the wings 
were swept back by shearing the cross sections parallel to 
the air stream, the span and aspect ratio remained constant. 
For the other case, a series of rectangular-plan-form wings 
were mounted on a special base which could be rotated to 
provide any desired angle of sweepback. This rotatory 
base was also used to examine the critical speed of swept- 
forward wings. 

Tests were conducted also on special models that were of 
the “ rotated ’' 1 ty^pe (sections normal to the leading edge 
were the same at all sweep angles) with the difference that 
the bases were alined parallel to the air stream. Two series 
of such rotated models having different lengths were tested. 

Inasmuch as the location of the center of gravity, the mass- 
density ratio, and the Mach number have important effects 
on the flutter characteristics of unswept wings, these param- 
eters were varied for swept wings. In order to investigate 
possible changes in flutter characteristics which might be 
due to different flow over the tips, various tip shapes were 
included in the experiments. 

In an analysis of flutter, vibrational characteristics are 
very significant; accordingly, vibration tests were made on 
each model. A special study of the change in frequency 
and mode shape with angle of sweep was made for a simple 
aluminum-alloy beam and is reported in appendix A. 

Theoretical analyses of the effect of sweep on flutter exist 
only r in brief or preliminary forms. In England in 1942, 
W. J. Duncan estimated by certain dimensional consider- 
ations the effect of sweep on the flutter speed of certain 
specialized iving ty-pes. Among other British workers 
whose names are mentioned in connection with problems 


J Supersedes NACA TN 2121, “Study of Effects of Sweep on the Flutter of Cantilever Wings” by J. G. Barmby, H. J. Cunningham, and I. E. Garrick, 1950. 
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of flutter involving sweep are R. McKinnon Wood, A. R. 
Collar, and I. T. Minhinnick. An account of Minhinnick’s 
work was given by Broadbent in reference 2. In reference 3 
a preliminary analysis for the flutter of swept wings in 
incompressible flow is developed on the basis of a “ strip 
theory’’ (with the strips taken in the stream direction) and 
is applied to the experimental results of reference 1. Exam- 
ination of the limiting case of infinite span discloses that 
the aerodynamic assumptions employed in reference 3 
are not well-grounded. Reference 4 adapts this strip 
theory to flexible wings and also presents an alternative 
“velocit} r component” treatment employing other aerody- 
namic assumptions which in their end result appear more 
akin to those employed in the analysis of the present report. 
No definite choice is made in reference 4 between the two 
methods although the strip-theory method is favored. 

In the present report a theoretical analysis is developed 
anew and given a general presentation. Application of the 
analysis has been limited at this time chiefly to those calcu- 
lations needed for comparison with experimental results. 
A wider examination of the effect of various parameters 
and of additional degrees of freedom on the flutter character- 
istics is desirable. 

SYMBOLS 


6 

b r 

V 


c 

l 

A 

A 

x' 

y' 

z' 

z 

v 

£ 

h 

e 


a 


My'), FhM 
My'), 

t 

CO 


half-chord of wing measured perpendicular to 
elastic axis, feet 

half-chord perpendicular to clastic axis at 
reference station, feet 

effective length of wing, measured along 


elastic axis, feet- 

wing chord measured perpendicular to elastic 
axis, inches 

length of wing measured along midchord line, 
inches 

angle of sweep, positive for sweepback, degrees 

geometric aspect ratio 

coordinate perpendicular to elastic axis in 
plane of wing, feet 
coordinate along elastic axis, feet 
coordinate in direction perpendicular to 
x'y'- plane, feet 

coordinate of wing surface in ^'-direction, feet 
nondimensional coordinate along elastic axis 

(y'ln 


( (l cos A) 2 ^ 


coordinate in wind-stream direction 
bending deflection of elastic axis, positive 
downward, feet 

torsional deflection of elastic axis, positive 
with leading edge up, radians 

local bending slope of elastic axis (^) 


local rate of change of twist 


W/ 


deflection function of wing in bending 
deflection function of wing in torsion 


time 

angular frequency of vibration, radians per 
second 


0>h 


co a 


A 

A 

ft 

fa 

fe 

f R 

f A 

V 

Ve 

Vn 

V e 

Vr 

Vr 

Fa 

V D 

k n 


<P 


P 


<Z 

M 

M e 

X cs 


X ea 



a 

a-\-x a 


m 

IC 

la 


angular uncoupled bending frequency, radians 
per second 

angular uncoupled torsional frequency about 
elastic axis, radians per second 

first bending natural frequency, cycles per 
second 

second bending natural frequency, cycles 
per second 

first torsion natural frequency, cycles per 
second 

uncoupled first torsion frequency relative to 
elastic axis, cycles per second 



(xjr a y y\ 

i-ca/a> 2 J ; 


experimental flutter frequency, cycles per 
second 

reference flutter frequency, cycles per second 
flutter frequency determined by analysis of 
present report, cycles per second 
free-stream velocity, feet per second 
experimental flutter speed, feet per second 
component of air-stream velocity perpen- 
dicular to elastic axis, feet per second 
( v cos A) 

experimental flutter speed taken parallel to 
air stream, miles per hour 
reference flutter speed, miles per hour 
reference flutter speed based on wing elastic 
axis, miles per hour (defined in appendix B) 
flutter speed determined by theory of present 
report, miles per hour 

theoretical divergence speed, miles per hour 
reduced frequency employing velocity com- 
ponent perpendicular to elastic axis (c cb/v n ) 
phase difference between wing bending and 
wing torsion strains, degrees 
density of testing medium at flutter, slugs per 
cubic foot 

dynamic pressure at flutter, pounds per square 
foot 

Mach number at flutter 
critical Mach number 

distance of center of gravity behind leading 
edge taken perpendicular to elastic axis, 
percent chord 

distance of elastic center of wing cross section 
behind leading edge taken perpendicular to 
elastic axis, percent chord 
distance of elastic axis of wing behind leading 
edge taken perpendicular to elastic axis, 
percent chord 


nondimensional elastic-axis position 



o 


nondimensional center-of-gravity position 



mass of wing per unit length, slugs per foot 
wing mass-density ratio at flutter (tt pb 2 /m) 
mass moment of inertia of wing per unit length 
about elastic axis, slug-feet“ per foot 


STUDY OF EFFECTS OF SWEEP ON THE FLUTTER OF CANTILEVER WINGS 


3 


r a 


nondimensional radius of gyration of wing 


about clastic axis 



El bending stiffness, pound-inches 2 in tables, 

pound-feet 2 in analysis 

GJ torsional stiffness, pound-inches 2 in tables, 

pound-feet 2 in analysis 

g h structural damping coefficient for bending 

vibration 

g a structural damping coefficient for torsional 

vibration 

P oscillatory lift per unit length, positive down- 

ward (defined in equation (6)) 

Ma oscillatory moment about elastic axis, positive 

leading edge up (defined in equation (7)) 

( } a special bracket used to identify terms which 

are due solely to inclusion of the last term 
in equation (5b) 

In order to preserve continuity and to facilitate comparison 
with previous work on the unswept wing, the subscript a 
rather than 6 is retained with certain quantities to refer to 
the torsional degree of freedom. 


ANALYTICAL INVESTIGATION 

GENERAL 

Assumptions. — An attempt is first made to point out the 
main assumptions which seem to be applicable for swept 
wings of moderate taper and of high or moderate length- 
chord ratios. 

(a) The assumptions, such as small disturbances and poten- 
tial flow, commonly employed in linearized treatment of 
unswept wings in an ideal incompressible fluid are made. 

(b) The structural behavior is such that over the main part 
of the wing the elastic axis may be considered straight . The 
wing is also considered sufficiently stiff at the root so that it 
behaves as if it were clamped normal to the elastic axis. 
An effective length V needed for integration reasons may be 
defined (for example, as in fig. 1). The angle of sweepback 
is measured in the plane of the wing from the direction nor- 
mal to the air stream to the elastic axis. All section param- 
eters such as semichord, locations of elastic axis and center 
of gravity, radius of gyration, and so forth, are based on 
sections normal to the elastic axis. 



Figure 1 .— Nonuniform swept wing treated in the present analysis. 


(c) The aerodynamic behavior is such that any section dy' 
of the wing normal to the elastic axis, taken in the direction 
of the component v cos A of the main-stream velocity, gener- 
ates a velocity potential associated with a uniform infinite 
swept wing having the same instantaneous distribution over 
the chord of velocity normal to the wing surface as does the 
actual section. 

Additional remarks on these assumptions are appropriate. 
With regard to assumption (a), in accordance with lineariza- 
tion of the problem, the boundary conditions are stated and 
treated with respect to a reference surface, in this case a 
plane, containing the mean equilibrium position of the wing 
and the main-stream velocity. Furthermore, incompressible 
flow is assumed in order to avoid complexity of the analysis, 
although modifications due to Mach number effects can be 
added. Such modifications may be based, for example, for 
wings having large length-chord ratios, on existing theoretical 
calculations of aerodynamic coefficients for subsonic or super- 
sonic two-dimensional flow appropriate to the component 
v cos A. On the other hand the modifications may be partly 
empirical, especially for “transonic” conditions and for small 
length-chord ratios. The transonic conditions and the gen- 
eral aerodynamic behavior of swept wings may depend, for 
large length-chord ratios, on the component v cos A, but the 
dependence may shift to the stream velocity v for small length- 
chord ratios. 

With respect to assumption (b), results of analyses of and 
experiment on unswept wings having low ratios of bending 
frequency to torsion frequency show that small variations of 
position of the elastic axis are not important. The assump- 
tion of a straight elastic axis over the main part of a swept 
wing, similarly, is not critical for many cases. This assump- 
tion is made for convenience, however, and modifications for 
a curved elastic axis can be made when necessary, for example, 
for plate-like wings. Small differences in the angle of sweep- 
back of the leading edge, quarter-chord line, elastic axis, and 
so forth, are neglected. The analysis could be further modi- 
fied to take into account variation of the angle of sweepback 
along the length of the wing. 

Assumption (c) implies that associated with the action 
of the wing in pushing air downward there is a noncirculatory 
potential-type flow similar to that around sections of an 
infinite flat-plate wing. Furthermore, as in the case of the 
unswept airfoil, a circulatory potential-type flow is generated 
in which for the swept airfoil the component v cos A is 
decisive in fixing the circulation. (This assumption differs 
from that made in the strip theory of references 3 and 4 
which employs the main-stream velocity together with 
sections of the wings parallel to the stream direction.) 
Effects of the floating of the wake in the stream direction 
rather than in the direction of v cos A and induced effects 
of variation of the strength of the wake in the wing-length 
direction are neglected, as are three-dimensional tip effects. 
For large values of the reduced frequency k n , a given segment 
of the wing might be influenced chiefly by the nearby wake 
and the correction would be small. On the other hand, for 
small values of k n a given segment might be influenced by a 
more widespread portion of the wake; corrections for this 
condition may possibly be based on knowledge of the static 
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case (for example, slope of the lift curve). As the angle of 
sweep approaches 90°, obviously the mechanism for the 
generation of lift is different from the one postulated here; 
for example, a tip condition may replace the trailing-edge 
condition and considerations of very small aspect ratio arise. 

Basic considerations. — Consider the configuration shown 
in figure 1 where the vertical coordinate of the wing surface 
is denoted by z' — Z(x' ,y' ,t) (positive downward). The 
effect of the position and motion of the wing may be given 
by the disturbance-velocity distribution to be superposed 
on the uniform stream in order to represent the condition 
of tangential flow at the wing surface. This velocity 
distribution normal to the surface (positive upward) is, for 
small disturbances, 


involving 0 and r vary in a linear manner. The noneircula- 
tory velocity potentials as in reference 5 and the new poten- 
tials associated with a and r are 


<t>h = hb yj 1 — X 2 

<t>o=v n Ob yl — x 2 
<t>o=v n (i tan \b\\— x 2 

4> e — Qb 2 (^ — a) V'l -T 2 
<t> r =v n T tan A b 2 (j — a) ^l—x 2 


(4) 


w(x',y'.t)=^ t +v-~ (la) 

where £ is the coordinate in the wind-stream direction. 
With the use of the relation 


dZ dZ dx' dZ dy' 
d£~d? d£ ' by' d£ 

dZ az . 

= x , cos A + . , sin A 
dx dy 

the vertical velocity at any point is 

, , , .v dZ dZ dZ. , , 

w(x',y . t) =~^ + v ^ cos A + v , sin A (1 b) 

Let the wing be bending so that a segment dy' (see fig. 1) is 
displaced from its equilibrium position by an incremental 
distance h (positive down) and also let the wing segment be 
twisting about the elastic axis through an incremental angle 0 
(positive leading edge up). The position of each point of 
the segment may be defined, for small deflections, by 

Z—h + x'O (2) 

The velocity distribution normal to the surface, equation (lb), 
conseq uen t ly l )ecomes 

W=h-\~: r'0 + vO cos A+ v(a + x'r) sin A (3) 


where cr=^, is the local bending slope of the elastic axis 

d9 

and is thus analogous to dihedral, and where T= ^ } y 


local change of twist of the elastic axis. 

In accordance with assumption (c) the noncirculatory-flow 
velocity potentials associated with the vertical- velocity dis- 
tribution are first needed. In equation (3) the terms involv- 
ing h , 0 , and a are constant across the chord, whereas those 


where v n = v cos A and x is the nondimensional chordwise 
coordinate measured from the midchord as in reference 5 
and related to the coordinate x' in the manner 

x =J+a 

The velocity potential for the circulatory flow associated 
with the wake may be developed on the basis of assumption (c) 
and the concepts for the infinite unswept wing introduced 
in reference 5. (Thus the circulatory-flow pattern for a sec- 
tion dy' of the finite swept wing is to be obtained from the 
corresponding flow pattern for an infinite uniform yawed 
wing. This infinite wing is assumed to have undergone har- 
monic oscillations for a long time; the full wake is established, 
remains where formed, and consequently is harmonically 
distributed in space. For the infinite uniform yawed wing, 
results for the circulatory flow are like those of reference 5 
with v replaced by the component v n and with the addition 
of terms to take care of o and r.) In particular, the strength 
of the wake acting on each section is determined by the condi- 
tion of smooth flow (the velocity remaining finite) at the trail- 
ing edge. This condition is utilized in the form ^ (<f>r + 0tf) 

is equal to a finite quantity at the trailing edge (where </>r is 
the velocity potential due to the vorticity in the wake, and 
(f> N is the total noncirculatory velocity potential), and this 
condition leads to a relation analogous to equation (VII) of 
reference 5 involving the basic quantity 

Q=h-\-v n 0+v n a tan A +6 (^- — o}j(0-\-v n T tan A) 

which occurs in the terms associated with the wake. The net 
result of these considerations is that the circulatory-flow 
velocity potential may be regarded as determined. 


STUDY OF EFFECTS OF SWEEP ON THE FLUTTER OF CANTILEVER WINGS 


5 


The pressure difference between upper and lower surfaces of the wing at a point x is (positive downward) 


p= — 2 P 


= — Zp 



d</A 


. , d(f>. 

v ^7 cos A+ V ^7 sin A 


) 


(5a) 


where 0 is in general the total potential (the sum of circulatory-flow and noncirculatory-flow potentials). The last term in 
equation (5a) is the product of the component of main-stream velocity taken along the wing and the lengthwise change 
in the velocity potential and is often neglected even in steady-flow work. The question of the retention or neglect of this 
last term seems partly dependent on the order in which the approximations are introduced — specifically, whether velocity 
potentials for the whole flow pattern are found and then the integrated forces are determined or whether section forces are 
first determined and then integrated. It seems appropriate to retain at least the noncirculatory part <t> N of 0 in the last 
term of equation (5a). In view, however, of the nature of the approximate treatment of the circulatory potential and of the 
inherent shortcomings of a strip analysis, in particular the neglect of lengthwise variations in wake vortex strength, compli- 
cating the results by also including 0 r in this term does not appear worth while. (This neglect of 0 r and retention of cj) N is 
realized to involve some inconsistencies in that account may not be taken of other higher order terms associated with length- 
wise variation of the wing wake, which may be of the same order as terms retained.) Thus equation (5a) becomes 


V~~ 


-2 P 


(*± +v 

V dt ' dx 


) , C )<t>N . \ 

7 cos A + sin A J 


(5b) 


For harmonic motion in each degree of freedom, relations for the pressure may be integrated over the chord to yield 
expressions for the air forces and moments. For the sake of separating and identifying the terms in force and moment ex- 
pressions which are due solely to the inclusion of the last term in equation (5b), a special bracket { } is employed. Thus 

these terms may be readily omitted. Numerical checks among the calculations made for the present report showed the 
effect of inclusion of the last term in equation (5b) on the calculated results to be quite small, even for 60° of sweepback 
within the range of other parameters investigated. 

The expressions for the aerodynamic lift (positive down) and for the moment about the elastic axis (positive leading 
edge up), each per unit length of the wing, are as follows: 


P= — 2irpv n bC ^Ji + Vnd+Vna tan A +b ( 0 +r n r tan A)J— 

Tpb^h + v n d + v n a tan A+[v n cr tan A + v n 2 r tan A + v n 2 tan 2 AjJ+ 

7 r pb 3 a + tan A+^y rt f tan A-\-v n 2 tan 2 aIJ 

M a = 2Tpv n b 2 O + + tan A + 6 ( d-\-v n T tan A)J— 

TTpv n b 3 0+^ v n T tan aJ+ttp/; 3 c/ j^A+ y n o- tan A + 

yv n (j tan A + vJt tan A + v n 2 tan 2 A^J— 7rp6 4 j^0 + 

A-f^ n f tan A-\-v n 2 tan 2 A^J 


(6) 


v n T tan A 


(7) 


6 

where 
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C =C(k n )=F(k n ) + iG(k n ) 


is the function associated with the wake developed by 
Theodorsen in reference 5; the reduced frequency parameter 
k n is defined by 


k n 


c ob 

V n 


cob 

v cos A 


(8) 


As has already been stated, the foregoing expressions were 
developed and apply for steady sinusoidal oscillations, 


6 ~ di (y f ) 6 iu>t 


(9) 


The amplitude, velocity, and acceleration in each degree of 
freedom are related as in the degree of freedom h\ that is, 


h = icoh 
h = — co 2 h 


Expressions for force and moment. — With the use of such 
relations, equations (6) and (7) may be put into the form 


P = — 7rp b 3 co 2 (h B ch -f- BB c e) 
M a = — 7rpb A co 2 (h B ah + OB a o) 


(10a) 

(11a) 


where 


B ch =^A ch +^ tan A (-i A) ({-1J + A0 + 
B„=A„+^ b tan A * (~p)) 


B a fi — i A a h~\~ 


tan A (— i ^ ) ((a) + /!„),) 


+ 


{x w A ( _ i6)} 

* » <»,' a w->+ A (!+»•) &■) 


FqQ — (i a I 


in which the four following coefficients: 

2G..2F 

A C h= — i— j — hi -j— 

1C n 1C n 


A.=a+p-(|-«) ff+i [e+JEP+G -0 ) I!] 
a,-«+G+«) “+« (f +«) (-0 

= -^-Q+a) A ch 

A*=- a 2 -(i+a) f5+(i- a2 ) U + 


are identical with those used in the case of the unswept wing. 
Additionally, 


A,= “^[ _ 6 + {^ +a2 }) +i ^G _lf ' 1 ) _ G _a2 )A A ] 


It is of interest to note that equations (6) and (7) reduce, 
for the case of the wing in steady flow (k n = 0), to 

P= — 2Tpbv n 2 1^0+ <r tan A+ t b tan A + 9 — + 

{2 w taa2 A ~l 62 w tan2 A }] (10b) 

M a = 2Tpb 2 v n 2 [(.+„ tan A)(l+a)+ra6 tan a(Q^J — a)+ 

{t w tan2 62 G +a2 ) w tan2 A }] (llb) 


per unit length of wing. 

Introduction of modes. — Equations (10a) and (11a) give 
the total aerodynamic force and moment on a segment of a 
sweptback wing oscillating in a simple harmonic manner. 
Relations for mechanical equilibrium applicable to a wing 
segment may be set up, but it is preferable to bring in directly 
the three-dimensional-mode considerations. (See for example, 
reference 6.) This end may be readily accomplished by 
the combined use of Rayleigh type approximations and the 
classical methods of Lagrange. The vibrations at flutter 
are assumed to consist of a combination of fixed mode shapes, 
each mode shape representing a degree of freedom associated 
with a generalized coordinate. The total mechanical 
energy, the potential energy, and the work done by applied 
forces, aerodynamic and structural, are then obtained by the 
integration of the section characteristics over the span. 
The Rayleigh type approximation enters in the representa- 
tion of the potential energy in terms of the uncoupled 
frequencies. 

As is customary, the modes are introduced into the problem 
as varying sinusoidally with time. For the purpose of sim- 
plicity of analysis, one bending degree of freedom and one 
torsion degree of freedom are carried through in the present 
development. Actually, any number of degrees of freedom 
may be added if desired, exactly as with an unswept wing. 
Let the mode shapes be represented b}^ 


h = [W)\k | 

d=[fe(y')\d | 


(12) 


where h=h 0 e i(al is the generalized coordinate in the bending 
degree of freedom, and 9_=d 0 e io}t is the generalized coordinate 
in the torsion degree of freedom. (In a more general treat- 
ment the mode shapes must be solved, but in this procedure 
f h (y') and jo{y') are chosen, ordinarily as real functions of y'. 
Complex functions may be used to represent twisted modes.) 
The constants h„ and d a are in general complex and thus 
signify the phase difference between the two degrees of 
freedom. 
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In the subsequent treatment the reader will notice that 
in some expressions, namely for force and moment, h and 8 
can conveniently and logically be retained in their complex 
form. In other expressions, notably for energy, one is 
forced to utilize h and 6 as real quantities. Appropriate 
statements will be made where necessary. 

For each degree of freedom an equation of equilibrium 
may be obtained from Lagrange's equation 


^/dT\_dT dC7 = 
dt\dq { ) dqidqt * l 


(13) 


where q t is a generalized coordinate and Q t is the correspond- 
ing generalized force. The kinetic energy of the mechanical 
system is 


T=l + l 6 2 f‘h {je(y')W + 

h,6_j o mx a b{j h (y')\ [je(y')]dy' (14) 


modes the frequencies appropriate to pure modes (obtained 
by proper constraints) are often used. On the other hand, 
employment of the normal or natural modes and frequencies 
appropriate to them, which might be obtained by proper 
ground test or by calculation, may be preferred. In either 
case the convenience of not having cross-stiffness terms in 
the potential-energy expression is noted. 

Application is now made to obtain the equation of equilib- 
rium in the bending degree of freedom. Equation (13) 
becomes 


d 

Jt 




dU 

dh 


— Qh 


( 16 ) 


The term Q h represents all the bending forces not derivable 
from the potential -energy function and consists of the aero- 
dynamic forces together with the structural damping 
forces. The virtual work 8W done on the wing by these 
forces as the wing moves through the virtual displacements 
5 h and 86 is 


where h and 0 here and in the subsequent equations (15) are 
to be interpreted as real in order that the energy be always 
positive (or zero), and for definiteness can be regarded as the 
real parts of h 0 e iut and 6 o e i03t respectively; and where 
m mass of wing per unit length, slugs per foot 
I a mass moment of inertia of wing about its elastic axis 
per unit length, slug-feet 2 per foot 
x a b distance of sectional center of gravity from the elastic 
axis, positive rearward, feet 

The potential energy of the mechanical system may be 
expressed in a form not involving bending-torsion cross- 
stiffness terms: 



' v+ kC ej (w) ,v 


(15a) 


where 

El bending stiffness, pound-feet 2 

GJ torsional stiffness, pound-feet 2 

If Rayleigh type approximations are used to introduce 
frequency, the expression for the potential energy may be 
written in a more convenient form: 


U =•! <*hk 2 ^ mf n 2 dy'+^ w« 2 0 2 J^ I«je 2 dy' (15b) 
Another expression for the potential energy is 

u=\ A 2 JJ C h j,Hy'+\ i f‘' CJ» 2 dy' (15c) 

The effective spring constants C h and C a correspond to unit 
length of wing and thus conform to their use in references 
5 to 7. The constants are effectively defined bv 



5 W = J o ‘ ' [(/>- C h ^ A) 8h+(M a - C a £ *)««] dy’ 

— f (P—mu h 2 ^J h h^fndy'8h+ 

£ (M a - I a u a 2<J ^M).fody' 8d=Q h 8h+ Q e 86 (17) 

where 

g h structural damping coefficient for bending vibration 

g a structural damping coefficient for torsional vibration 

In this expression the aerodynamic forces appropriate to 
sinusoidal oscillations are used. The application of the 
structural damping as in equation (17) (proportional to 
deflection and in phase with velocity) corresponds to the 
manner in which it is introduced in reference 7. In accord- 
ance with the preceding development, the aerodynamic and 
structural damping forces and moments in equation (17) 
are regarded as complex, but the virtual displacements 8h 
and 88 should be considered real. Thus, the physically 
significant part of the resulting expression for virtual work is 
the real part. Since the subsequent analysis reverts to ex- 
pressions for forces and moments, no further qualifications on 
the use of h and 6 in their convenient complex forms are 
needed. 

F or the half-wing 
Qh= (P- m^^U^f.dy' 

^ cft ) (tan A)/,^';+ JA^ b (tan 2 A)/ A 0d + 

6A ca je f h -\-dA cr b (tan A)f h + 


These effective spring constants are related to the frequencies 
associated with the chosen modes. For so-called uncoupled 

042523-52 2 



(^) <lijif\dy' (18) 
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where b r is the semichord at some reference section. Performance of the operations indicated in equation (16) and collec- 
tion of terms lead to the equation of equilibrium in the bending degree of freedom 



where 


>Vl 1 r r /b \ 2 t' 1 ' / bV / 1 \ 

-) - .fndy' — j- ( ,) A cl J h 2 dy'+ (t-) tan A( i )({— 1} + 

r/ K UrJ 0 \0t/ JO \0 r / \ ^n/ 


rr / b x 

b r \ ( t~ ) Act (tan 


(19a) 


1 m 

K Tpb~ 


By a parallel development the equation of equilibrium for the torsional degree of freedom may also be obtained as follows: 

(*{c f '(£)‘( t- - 4 -*) ^+X'(t)‘*“ A 0£) {M+A - )U % dy,+ { "'iXxf A (£?) *• P d A\ 

•- {[' -(5)‘«+«#4 JT(0 T w-fXr)’ *">***- 6 'f (0 tan A( ' w/ * t *<- 
{‘■'X* (£)* tan ' A (s + “’) £■ '• ' w J I) 


+ 


(20a) 


where r a = V / Jmb 2 (radius of gyration of wing about the elastic axis). 

Determinantal equation for flutter. — Equations (19a) and (20a) may lie rewritten with the use of the nondimensional 


coordinate 7 l=jr' They then are in the form 


where 


(hA 2 -]-&^2) 7r P bjcx) 2 = 0 

2~\~ 2) 7T p b r A 0}~— 0 


(19b) 

(20b) 


X ’ 0 a (i a- 1 ) + ^ [f.< rfr, - {^X' t°^A( 

c io//)\ 3 /r \ f 1H //>\ 4 <IFb 

'*=/' J o ( yj (• A ea ) [F„(y)} [F,(„)l dy - b r j (A) tan A (A CT ) [F h ( ,)] dy + 

tbr f 10 / 6 V, 21 / a \ r t’ / m ^ 2 f<> , 1 

t /' Jo (.6 J 1 A \l-„ 2 r dy 2 ^ 1 


//, 


y ri.o/ /)\3 /r \ ri.o// ) \-t / 1 \ ///<■. 

= [ j o (£ ) ( ; - A*) [F,(„)] [F„(y)\ dy + J ( A ) tan A (i j- ) ({a] + A ah ) [F»(y)] }[ A dy + 


dF h 


[Fe(y)\ <l £* dy 


} 


r \ 2 n /*i.o/ 7 ) \4 v. 2 /m .0/ \ 4 

(A) ^-[F^)] 2 ^-rJ o (■^~) A a a[F«{y)] 2 dy 

IFo , tbr 2 f 10 / b \® 2 /I . A 1 d a F, , 1 

?v J. w A (>+*•.) f? inwi ./,■ *] 


/ X tan A(d ar ) [F,(i,)] ~ 


in which F h (rj) = f h (l f r}) and F e {r])=je{l f t)). 

The borderline condition of flutter, separating damped and undamped oscillations, is determined from the nontrivial 
solution of the simultaneous homogeneous equations (19b) and (20b). Such a solution corresponds to the fact that median- 
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ical equilibrium exists for sinusoidal oscillations at a certain airspeed and with a certain frequency. The flutter condition 
thus is given by the vanishing of the determinant of the coefficients 


A 2 

D 2 


b 2 


= 0 


E s 


Application to the case of uniform, cantilever swept wings is made in the next section. 


APPLICATION TO UNIFORM CANTILEVER SWEPT WINGS 

The first step in the application of the theory is to assume or develop the deflection functions to be used. For the purpose 
of applying the analysis to the wing models empIo} r ed in the experiments it appeared reasonable to use for the deflection 
functions, F h (rj) and F e (rj), the uncoupled first bending and first torsion mode shapes of an ideal uniform cantilever beam. 
Although approximations for these mode shapes could be used, the analysis utilized the exact expressions developed from 
equations (120) and (106d), respectively, of reference 8 by application of appropriate boundary conditions. 

The bending-mode shape can be written 


F ^= C ' [ rsMl+cosgl (cos/^-cosh/M+sinh 


-sin 


■’] 


where 0i = O.59697r for first bending. The torsion mode shape can be written 

F e (r))=C 2 sin i Q 2 r] 

where /3 2 =-a for first torsion and C x and C 2 are constants. 

The integrals appearing in the determinant elements A 2 , B 2 , D 2) and E 2 are 


i.° (IF, 

F h dr ) = 3.7110C l 2 

1,0 F t <^ = 0.3183(7 2 2 


I F h 2 dij = 1.8554 Ci 2 

Jo 

Fe 2 dti = 0.5000C2 2 

f 

x: 
x 

f F h Fedni = f F e F h dr) = — 0.$233C 1 C 2 
Jo Jo 

X 

x; 


.o d 2 F k 

F„ a ^»d v = 1.5926C l * 


i.o ,n F a 

F e J rfjj = -1.2337C , 2 2 


i.o d Fa 

F h dt) — — \ .4040C'iC 2 


L0 r ]F, 

F,~ d v = -2.0Q&9C 1 C 2 


ri.o ( pFa 

J () F h (l J dr, = 2.27 82C,C 2 
i.o ( PF, 

F e dt\ = — \A122G\C 2 


f 


The flutter determinant becomes 


1.8554C,* j- -4 + 3.7110C,* (* £-) ([-lJ + .-l.*) tan A.- ^1.592 6C, 2 
-0.92330,0^- W-2.0669 C i C 2 ^'-) ([a) + .l„,,) tan A+ 1.4722C,C S 
or more conveniently, when columns and rows of the determinant are divided by appropriate terms 


— 0.9233C, C 4 t'ii - ( - 1 .4040C, C,) b r A c , tan A + 1 2.2782C, Cjb, ~ 
0.o000C,*l'f?-0.3183C,*6, .1., tan A-f- l.2337C a J i> r ' (|+° 2 )/ a } 


= 0 




0.9189 D +2.0571 1 ‘" 1 


where 


wtOk) +{ , 4 * 5 *( , w)' &) E - o ' 63660 ^^+( 2 ' 467 <^;') , (l +< ‘ , )i?} 


=0 


B=~— A c 


D=——A ah 


e = t [ 1 -(^) 2(l +*?«)]- ^ 


It is interesting to note that the parameters A and V jb T appear only in the combination * n the immediately preceding determinant . The solution of 

the determinant results in the flutter condition. 
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EXPERIMENTAL INVESTIGATION 
APPARATUS 

Wind tunnel. — The tests were conducted in the Langley 
4.5-foot flutter research tunnel which is of the closed-throat, 
single-return type employing either air or Freon-12 as a 
testing medium at pressures varying from 4 inches of mercury 
to 30 inches of mercury. In Freon-12, the speed of sound 
is 324 miles per hour and the density is 0.0106 slug per cubic 
foot at standard pressure and temperature. The maximum 
choking Mach number for these tests was approximately 
0.92. The Reynolds number range was from 0.26X10 6 to 
2.6X10 6 with most of the tests at Reynolds numbers of the 
order of 1.0 X10 6 . 

Models. — In order to obtain structural parameters re- 
quired for the flutter studies, different types of construction 
were used for the models. Some models were solid spruce, 
others were solid balsa, and many were combinations of 
balsa with various aluminum-alloy inserts. Seven series of 
models were investigated, for which the cross sections and 
plan forms are shown in figure 2. 

Figure 2 (a) shows the series of models which were swept 
back by shearing the cross sections parallel to the air stream. 
In order to obtain flutter with these low-aspect-ratio models, 
thin sections and relatively light and weak wood construction 
were employed. 



''Lam/noted spruce 
8 " 7 s - 



Models 11-15 

Sections parallel to the air stream 
(a) Sheared swept models with a constant geometric aspect ratio of 2. Series I. 
Figure 2. — Model plan form and cross-sectional construction. 



(D) Models swept back by use of a rotating mount. Series II. 
Figure 2.— Continued. 


The series of rectangular-plan-form models shown in figure 
2 (b) were swept back by using a base mount that could be 
rotated to give the desired sweep angle. The same base 
mount was used for testing models at forward sweep angles. 
It is known that for forward sweep angles divergence is 
critical. In an attempt to separate the divergence and 
flutter speeds in the sweepforward tests, a D-spar cross- 
sectional construction was used to get the elastic axis rela- 
tively far forward (fig. 2 (c)). 

Two series of wings (figs. 2 (d) and 2 (e)) were swept back 
with the length-chord ratios kept constant. In these series 
of models, the chord perpendicular to the leading edge 
was kept constant and the bases were alined parallel to the 
air stream. The wings of length-chord ratio 8.5 (fig. 2 (d)) 
were cut down to get the wings of length-chord ratio 6.5 
(fig. 2 (e)). 

Another series of models obtained by using this same 
manner of sweep (fig. 2 (f)) was used for investigating some 
effects of tip shape. 

Spanwise strips of lead were fastened to the models shown 
in figure 2 (e) and a series of tests were conducted with these 
weighted models to determine the effect of center-of-gravity 
shift on the flutter speed of swept wings. The method of 



Model 65 64 63 62 



(d) Swept models having a length-chord ratio of 8.5. Series IV 
Figure 2 — Continued. 
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varying the center of gravity is shown in figure 2 (g). In 
order to obtain data at zero sweep angle it was necessary, 
because ot the proximity of flutter speed to wing-divergence 
speed, to use three different wings. These zero-sweep-angle 
wings, of 8-inch chord and 48-inch length, had an internal 
weight system. 

The models were mounted from the top of the tunnel as 
cantilever beams with rigid bases (fig. 3). Near the root of 
each model two sets of strain gages were fastened, one set for 
recording principally bending deformations and the other 
set for recording principally torsional deflections. 


METHODS 

Determination of model parameters. —Pertinent geometric 
and structural properties of the model are given in tables I 
to VII. Some parameters of interest are discussed in the 
following paragraphs. 

As an indication of the nearness to sonic-flow conditions, 
the critical Mach number is listed. This Mach number is 
determined by the Karman-Tsien method for a wing section 
normal to the leading edge at zero lift. 

The geometric aspect ratio of a wing is here defined as 


a Semispan 2 _(l cos A) 2 _ l 2 A 
Plan-form area Ic c C ° S * 2 

The geometric aspect ratio A g is used in place of the con- 
ventional aspect ratio A because the models were only 
semispan wings. For sheared swept wings, obtained from a 
given unswept wing, the geometric aspect ratio is constant, 
whereas for the wings of constant length-chord ratio the 
geometric aspect ratio decreases with cos 2 A as the angle of 
sweep is increased. 

The weight, center-of-gravity position, and polar moment 
of inertia of the models were determined by usual means. 
The models were statically loaded at the tip to obtain the 
rigidities in torsion and bending GJ and EI. 

A parameter occurring in the methods of analysis of this 
report is the position of the elastic axis. A “section” elastic 
axis located at x ea was obtained for wings from each series 
of models as follows: The wings were clamped at the root 
normal to the leading edge and at a chosen spanwise station 


Mode/ 75 


74 


73 


72 


from trailing edge ot /" spacing 
Lengthwise do/so laminations. 




(e) Swept models having a length-chord ratio of 6.5. Series V. 
Figure 2.— Continued. 



(f) Models used to investigate the effect of tip shape on the flutter velocity. Series VI. 
Figure 2. — Continued. 



Model 

A 

(deg) 

91-1,91-2,91-3 * 

0 

92-1. 92-2. 92-3 

15 

93-1, 93-2, 93-3 

30 

94-1, 94-2, 94-3 

45 

95-1.95-2,95-3 

60 


* Chord =8", lead inside balsa. 

(g) Models used to determine the effect of center-of-gravity shift on 
the flutter velocity of swept wings. Series VII. * 

Figure 2.— Concluded. 



Figure 3.— Model 12 in the tunnel test section. 
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were loaded at points lying in the chordwise direction. The 
point for which pure bending deflection occurred, with no 
twist in the plane normal to the leading edge, was determined. 
The same procedure was used for those wings which were 
clamped at the root, not normal, but at an angle to the 
leading edge. A different elastic axis designated the “wing” 
elastic axis and located at x ea ' was thus determined. 

For these uniform, swept wings with fairly large length- 
chord ratios, the wing elastic axis was reasonably straight 
and remained essentially parallel to the section elastic axis, 
although it was found to move farther behind the section 
elastic axis as thi angle of sweep was increased. It is real- 
ized that in general for nonuniform wings — for example, 
wings with cut-outs or skewed clamping — a certain degree 
of cross stiffness exists and the concept of an elastic axis is 
an oversimplification. More general concepts such as those 
involving influence coefficients may be required. These 
more strict considerations, however, are not required here 
since the elastic-axis parameter is of fairly secondary impor- 
tance. 

The wing mass-density ratio k is the ratio of the mass of a 
cylinder of testing medium, of a diameter equal to the chord 
of the wing, to the mass of the wing, both taken for unit 
length along the wing. The density of the testing medium 
when flutter occurred was used in the evaluation of k. 

Determination of the reference flutter speed. -It is 
convenient in presenting and comparing data of swept and 
unswept wings to employ a certain reference flutter speed. 
This reference flutter speed will serve to reduce variations 
in flutter characteristics which arise from changes in the 
various model parameters such as density and section proper- 
ties not pertinent to the investigation. It thus aids in 
systematizing the data and emphasizing the desired effects 
of sweep including effects of aspect ratio and Mach number. 

This reference flutter speed V R may be obtained in the 
following way. Suppose the wing to be rotated about the 
intersection of the elastic axis with the root to a position of 
zero sweep. In this position the reference flutter speed is 
calculated by the method of reference 7, which assumes an 
idealized, uniform, infinite wing mounted on springs in an 
incompressible medium. For nonuniform wings, a reference 
section taken at a representative spanwise position, or some 
integrated value, may be used. Since the wings used were 
uniform, any reference section will serve. The reference 
flutter speed may thus be considered a “section” reference 
flutter speed and parameters of a section normal to the lead- 
ing edge are used in its calculation. This calculation also 
employs the uncoupled first bending and torsion frequencies 
of the wing (obtained from the measured frequencies) and 
the measured density of the testing medium at time of 
flutter. The calculation yields a corresponding reference 
flutter frequency which is useful in comparing the frequency 
data. For the sake of completeness a further discussion of 
the reference flutter speed is given in appendix B. 

Test procedure and records. — Since flutter is often a 
sudden and destructive phenomenon, coordinated test pro- 
cedures were required. During each test, the tunnel speed 
was slowly raised until a speed was reached for which the 


amplitudes of oscillation of the model in bending and torsion 
increased rapidly while the frequencies in bending and tor- 
sion, as observed on the screen of the recording oscillograph, 
merged to the same value. At this instant, the tunnel 
conditions were recorded and an oscillograph record of the 
model deflections was taken. The tunnel speed was im- 
mediately reduced in an effort to prevent destruction of the 
model. 

From the tunnel data, the experimental flutter speed V e , 
the density of the testing medium p, and the Mach number M 
were determined. No blocking or wake corrections to the 
measured tunnel velocity were applied. 

From the oscillogram the experimental flutter frequency 
f e and the phase difference <p (or the phase difference ±180°) 
between the bending and torsion deflections near the root 
were read. A reproduction of a typical oscillograph flutter 
record, which indicated the flutter to be a coupling of the 
wing bending and torsion degrees of freedom, is shown as 
figure 4. Since semispan wings mounted rigidly at the base 
were used, the flutter mode may be considered to correspond 
to the flutter of a complete wing having a very heavy fuselage 
at midspan — that is, to the symmetrical type. 

The natural frequencies of the models in bending and 
torsion at zero airspeed were recorded before and after each 
test in order to ascertain possible changes in structural char- 
acteristics. In most cases there were no appreciable changes 
in frequencies but there were some reductions in stiffnesses 
for models which had been weakened by fluttering violently. 
Analysis of the decay records of the natural frequencies 
indicated that the wing damping coefficients g h and g a 
(reference 7) were about 0.02 in the first bending mode and 
0.08 in the torsion mode. 


RESULTS ANI) DISCUSSION 


EXPERIMENTAL INVESTIGATION 

Presentation of experimental data. — Results of the 
experimental investigation are listed in detail in tables I to 
VII, and some significant experimental trends are illustrated 



Figure 4.— Oscillograph record of model at flutter. 
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in figures 5 to 10. As a basis for presenting and comparing 
the test results, the ratio of experimental tunnel stream 
conditions to the reference flutter conditions is employed so 
that the data indicate more clearly combined effects of aspect 
ratio, sweep, and Mach number. As previously mentioned, 
use of the reference flutter speed V R serves to reduce varia- 
tions in flutter characteristics which arise from changes in 
other parameters, such as density and section properties, 
which are not pertinent to this investigation. (See appen- 
dix B.) 

Some effects on flutter speed. — A typical plot showing 
the effect of compressibility on the flutter speed of wings 
at various angles of sweepback is shown in figure 5. These 
data are from tests of the rectangular-plan-form models 
(type 30) that were swept back by use of the rotating mount, 
for which arrangement the reference flutter speed does not 
vary with either Mach number or sweep angle. Observe 
the large increase in speed ratio at the high sweep angles. 

The data of reference 1 from tests of a rigid, flexibly 
mounted rectangular model having a rotating base are also 
plotted in figure 5. It can be seen that the data from the 
cantilever models of the present report which had a similar 
method of sweep are in conformity with the data from the 
flexibly mounted model. This indicates that, for uniform 
wings having the range of parameters involved in these 
tests, the differences due to mode shape are not very great. 

Figure 6 is a cross plot of the data from figure 5 plotted 
against A at a Mach number approximately equal to 0.65. 
The data of the swept wings of constant length-chord ratio 
and of the sheared swept wings are also included for comparison. 
The velocity ratio V e /V R is relatively constant at small sweep 
angles but rises noticeably at the large sweep angles. It is 
pointed out that the reference flutter speed V R may be 



for the rotated and constant-length-chord-ratio wings, but 
for the sheared wings this reference speed corres ponds to a 
curve decreasing somewhat less rapidly than Vcos A as A 
increases. (See appendix B.) 


The order of magnitude of some three-dimensional effects 
may be [noted from the fact that the shorter wings 

=6.5, fig. 6, series have higher velocity ratios than the 

longer wings ^=8.5, series iv) • This increase maybe due 

partly to differences in flutter modes as well as aerodynamic 
effects. 

Some effects on flutter frequency. —Figure 7 is a repre- 
sentative plot of the flutter-frequency data given in table II. 
The figure shows the variation in flutter-frequency ratio 
with Mach number for different values of sweep angle for 
the models rotated back on the special mount. The ordinate 


Series I n 3E 




0 10 20 30 40 50 60 70 

Sweep angle, A, deg 


Figure 6. — Cross plot of ratio of experimental to reference flutter speed as a function of 
sweep angle for various wings. Mach number is approximately 0.65. 



O ./ .3 .4 .5 .6 .7 .8 .9 

Mnrh number , , A/ 



0 ./ . 2 .3 .4 .5 .6 .7 .8 .9 

Mach number , A/ 


Figure 5.— Ratio of experimental to reference flutter speed as a function of Mach number 
for various sweep angles for series II models (fig. 2 (b)) on the rotating mount. 


Figure 7.— Ratio of experimental to reference flutter frequency as a function of Mach number 
for various sweep angles for series II models (fig. 2 (b)) on the rotating mount. 
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is the ratio of the experimental flutter frequency to the 
reference flutter frequency /*//*. It appears that there is a 
reduction in flutter frequency with increase in Mach number 
and also an increase in flutter frequency with increase in 
sweep. The data from reference 1 show the same trend 
with increase in sweep. Considerably more scatter may be 
noted in the frequency data than in the speed data (fig. 5) 
from the same tests. 

The results of the tests for rotated wings with cliordwise 
laminations (models 40A, B, C, D) are given in table II. 
At sweep angles up to 30° the values of the speed ratio 
V e JV R for wings of this construction were low (in the neigh- 
borhood of 0.9), and the flutter frequency ratios j e /j R were 
high (of the order of 1.4). As these results indicate and as 
visual observation showed, these models fluttered in a mode 
that apparently involved an appreciable proportion of the 
second bending mode. The models with spanwise lamina- 
tions (models 30A, B, C, D) also showed indications of this 
higher flutter mode at low sweep angles; however, these 
models were able to pass through the small speed range of 
higher mode flutter without sufficiently violent oscillations 
to cause failure. At a still higher speed these models with 
spanwise laminations fluttered in a lower mode resembling a 
coupling of the torsion and first bending modes. This 
lower mode type of flutter characterized the flutter of both 
the sheared and constant-length-chord-ratio models. 

For those wing models having the sheared type of balsa 
construction (models 22', 23. 24, and 25), the results are 
more difficult to compare with those of the other models. 
This difficulty arises chiefly because the lightness of the 
wood produced relatively high mass-density ratios k and 
partly because of the nonhomogeneity of the mixed wood 
construction. For high values of k the flutter-speed co- 
efficient changes rather abruptly even for the unswept models 
(reference 7). The data are nevertheless included in table I. 

Effect of shift in center-of-gravity position on the flutter 
speed of swept wings. — Results of the investigation of the 
effects of center-of-gravity shift on the flutter speed of 
swept wings are illustrated in figure 8. This figure is a 
cross plot of the experimental indicated air speeds as a 
function of sweep angle for various center-of-gravity posi- 
tions. The ordinate is the experimental indicated air 


speed 


r / P 
e \ 0.00238’ 


which 


serves to 


reduce the scatter 


resulting from flutter tests at different densities of testing 
medium. The data were taken in the Mach number range 
between 0.14 and 0.44, so that compressibility effects are 
presumably negligible. As in the case of unswept wings, 
forward movement of the center of gravity increases the 
flutter speed. Again, the flutter speed increases with 
increase in the angle of sweep. 

The models tested at zero sweep angle (models 91-1, 91-2, 
91-3) were of different construction from and of larger size 
than the models tested at the higher sweep angles. Because 
of the manner of plotting the results, namely as experimental 
indicated airspeed (fig. 8), a comparison of the results of 
tests at A=0° with the results of the tests of swept models 


is not particularly significant. The points at zero sweep 
angle are included, however, to show that the increase in 
flutter speed due to a shift in the center-of-gravity position 
for the swept models is of the same order of magnitude as 
for the unswept models. For the unswept models, the 
divergence speed V D and the reference flutter speed V R are 
fairly near each other, and although the models appeared 
to flutter, the proximity of the flutter speed to the divergence 
speed may have influenced the value of the critical speed. 

The method used to vary the center of gravity (see fig. 2 (g)) 
produced two bumps on the airfoil surface. At the low 
Mach numbers of these tests, however, the effect of this 
roughness on the flutter speed is considered negligible. For 
proper interpretation of figure 8, the fact must be kept in 
mind that the method of varying the location of the center 
of gravity changed the radius of gyration r a and the tor- 
sional frequency j a • 

The effect of sweepforward on the critical speed. — An 

attempt was made to determine the variation in flutter speed 
with angle of sweepforward by testing wings on the mount 
that could be rotated both backward and forward. As ex- 
pected, however, the model tended to diverge at forward 
sweep angles in spite of the relatively forward position of 
the elastic axis in this D-spar wing. 

Figure 9 shows a plot of the ratio of critical speed to the 
reference flutter speed V R against sweep angle A. Note the 
different curves for the sweptback and for the swept forward 
conditions and the sharp reduction in critical speed as the 


360 



40 


0 10 20 30 40 50 60 70 

Sweep angle , A, deg 

Figure 8.— Cross plot of flutter speed as a function of sweep angle for several center-of-gravity 
positions. Series VII models (fig. 2 (g) ). Length-chord ratio is approximately 6. 
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Figure 9. 


Sweep angle. A, deg 

-Comparison of sweepforward and sweepback tests on wings tested on a rotating 
mount. Series III models (fig. 2 (c)). 


angle of sweepforward is increased. The different curves 
result from two different phenomena. When the wing was 
swept back, it fluttered; whereas at forward sweep angles it 
diverged before the flutter speed was reached. Superposed 
on this plot for the negative values of sweep are the results 
of calculations based on an analytical study of divergence 
(reference 9). Reasonable agreement exists between theory 
and experiment at forward sweep angles. The small differ- 
ence between the theoretical and experimental results may 
perhaps be due to an inaccuracy in determining either the 
position of the elastic axis of the model or the required slope 
of the lift curve or both. 

The divergence speed V D for the wing at zero sweep angle, 
as calculated by the simplified theory of reference 7, is also 
plotted in figure 9. This calculation is based on the assump- 
tion of a two-dimensional unswept wing in an incompressible 
medium. The values of the uncoupled torsion frequency 
and the density of the testing medium at time of flutter or 
divergence are employed. Reference 9 shows that a rela- 
tively small amount of sweepback raises the divergence speed 
sharply. For convenience, however, the numerical quan- 
tity V D (based on the wing at zero sweep) is listed in table I 
for all the tests. 

Effect of tip modifications. — Tests to investigate some 
of the over-all effects of tip shape were conducted and some 
results are shown in figure 10. Two sweep angles and two 
length-chord ratios were used in the experiments conducted 
at two Mach numbers. It is seen that, of the three tip 
shapes used, namely, tips perpendicular to the air stream, 
perpendicular to the wing leading edge, and parallel to the 
air stream, the wings with tips parallel to the air stream gave 
the highest flutter speeds. 


DISCUSSION AND COMPARISON OF ANALYTICAL AND EXPERIMENTAL 

RESULTS 

Correlation of analytical and experimental results has been 
made for wings swept back in the two different manners; 
that is, (1) sheared back with a constant value of A e , and 
(2) rotated back. The two types of sheared wings (series I) 
and two rotated wings (models 30B and 30D) have been 
analyzed. 



Figure 10.— Effect of tip shape on the flutter speed of swept wings. Wings of length-chord 
ratios of 7.25 and 11 (fig. 2 (f)). Series VI models. 


Results of some solutions of the flutter determinant for a 
wing (model 30B) on a rotating base at several angles of 
sweepback are shown in figures 11 and 12. Figure 11 shows 
the flutter-speed coefficient as a function of the bending-to- 
torsion frequency ratio, and figure 12 shows the flutter fre- 
quency ratio as a function of the bending-to-torsion frequency 
ratio. 


The calculated results (for those wings investigated analyt- 
ically) are included in tables I and II. The ratios of experi- 
mental to analytical flutter speeds and flutter frequencies 
have been plotted against the angle of sweep in figures 13 
to 16. If an experimental value coincides with the corre- 
sponding analytical predicted value, the ratio will fall at a 
value of 1.0 on the figures. Deviations of experimental re- 
sults above or below the analytical results appear on the 
figures as ratios greater than or less than 1.0, respectively. 
The flutter-speed ratios plotted in figure 13 for the two ro- 
tated wings show very good agreement between analysis and 
experiment over the range of sweep angle, 0° to 60°. Such 
good agreement in both the trends and in the numerical 
quantities is gratifying but probably should not be expected 
in general. In view of the discussion of the last term in 
equation (5b) it may be of interest to mention that failure 
to include the terms arising from the last term of equation 
(5b) in the calculations for model 30B would decrease the 
ratio V e /V A corresponding to A=60° by about 3 percent. 
The flutter frequency ratios of figure 14 obtained from the 
same two rotated wings are in good agreement. 

The flutter-speed ratios plotted in figure 15 for the two 
types of sheared wings do not show such good conformit}^ at 
the low’ angles of sweep, whereas for sweep angles be}’ond 45° 
the ratios are considerably nearer to 1 .0. The sheared w ings 
are again observed to have a constant value of A g of 2.0 
(aspect ratio for the whole wing would be 4.0). For this 
small value of aspect ratio the finite-span correction is 
appreciable at zero angle of sweep and, if made, would bring 
better agreement at that point. Analysis of the corrections 
for finite-span effects on swept wings requires further 
consideration. 

Figures 13 and 15 also afford a comparison of the behavior 
of wings swept back in two manners: (1) rotated back with 
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Figure 11— Theoretical flutter-speed coefficient as a function of the ratio of bending to torsion 
frequency for the rotated model 30B at four angles of sweep and with a constant mass- 

density ratio ^-=37.8^ • 


constant length-chord ratio but decreasing aspect ratio (fig. 
13), and (2) sheared back with constant aspect ratio and 
increasing length-chord ratio (fig. 15). A study of these two 
figures suggests that the length-chord ratio rather than the 


aspect ratio 



may be the relevant parameter in 


determining corrections for finite swept wings. (Admit- 
tedly, effects of tip shape and root condition are also in- 
volved and have not been precisely separated.) 

Figure 16, which refers to the same sheared wings as figure 
15, shows the ratios of experimental to predicted flutter fre- 
quencies. The trend is for the ratio to decrease as the angle 
of sweep increases. Table 1 shows that the flutter frequency 
f R obtained with V R and used as a reference in a previous 
section of the report is not significantly different from the 
frequency /a predicted by the present analysis. 

A few remarks can be made on estimates of over-all trends 
of the flutter speed of swept wings. As a first consideration 
the conclusion may be made that, if a rigid infinite yawed 




Figure 13.— Ratio of experimental to theoretically predicted flutter speed as a function of 
sweep angle for two rotated models. 

wing were mounted on springs which permitted it to move 
vertically as a unit and to rotate about an elastic axis, the 
flutter speed would be proportional to 1/cos A. A finite 
yawed wing mounted on similar springs would be expected 
to have a flutter speed lying above the curve of 1/cos A 
because of finite-span effects. For a finite swept back wing 
clamped at its root, however, the greater degree of coupling 
between bending and torsion adversely affects the flutter 
speed so as to bring the speed below the curve of 1/cos A 
for an infinite wing. This statement is illustrated in figure 
17 which refers to a wing (model 30B) on a rotating base. 
The ordinate is the ratio of flutter speed at a given angle of 
sweep to the flutter speed calculated at zero angle of sweep. 
A theoretical curve is shown, together with experimentally 
determined points. Curves of 1/cos A and 1/Vcos A are 
shown for convenience of comparison. The curve for model 
30D (not shown in figure 17) also followed this trend quite 
closely. The foregoing remarks should prove useful for 
making estimates and discussing trends but are not intended 
to replace more complete calculation. In particular, men- 
tion may be made, for example, that a far-forward location 
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Fir. u re 14.- -Ratio of experimental to theoretically predicted (lutter frequency as a function 
of sweep angle for two rotated models. 



Figure 15. — Ratio of experimental to theoretically predicted flutter speed as a function cf 
sweep angle for two types of sheared models. 

of section center of gravity would lead to an entirely dif- 
ferent trend. Moreover, as is apparent from the analysis, 
the bending stiffness can play an increasingly significant 
role with increase in the angle of sweep. 

The experiments and calculations deal in general with 
wings having low ratios of natural first bending to first 
torsion frequencies. At high values of the ratio of bending 
frequency to torsion frequency, the position of the elastic 
axis becomes relatively more significant. Additional calcula- 
tions to develop the theoretical trends are desirable. 

CONCLUSIONS 

In a discussion and comparison of the results of an in- 
vestigation of the flutter of a group of swept wings, the 
manner of sweep is significant. This report deals with two 
main groups of uniform, swept wings: rotated wings and 
sheared wings. In presenting the data, employment of a 
certain reference flutter speed was found convenient. The 
following conclusions seem to apply: 

1. Comparison with experiment indicates that the analysis 
presented is satisfactory for giving the main effects of sweep, 
at least for nearly uniform cantilever wings of moderate 
length-chord ratios. Additional calculations are desirable 
to investigate various theoretical trends. 

2. The coupling between bending and torsion adversely 
affects the flutter speed. The fact, however, that only a 
part of the forward velocity is aerodynamically effective 



Sweep angle, A, deg 

Figure 17.— Flutter-speed ratio as a function of sweep angle for model 30B at a constant 
mass-density ratio showing analytical and experimental results. 

increases the flutter speed. Certain approximate relations 
can be used to estimate some of the trends. 

3. Although a precise separation of the effects of Mach 
number, aspect ratio, tip shape, and center-of-gravity posi- 
tion has not been accomplished, the order of magnitude of 
some of these combined effects has been experimentally 
determined. Experimental results indicated are 

(a) The location of the section center of gravity is an 
important parameter and produces effects for swept wings 
similar to those for unswept wings over the range (30 per- 
cent to 70 percent chord) of locations tested. 

(b) Appreciable differences in flutter speed have been 
found to be due to tip shape. 

(c) The length-chord ratio of swept wings is a more 
relevant finite-span parameter than is the aspect ratio. 

(d) Compressibility effects attributable to Mach number 
are fairly small, at least up to a Mach number of 0.8. 

(e) The sweptforward wings could not be made to flutter 
but diverged before the flutter speed was reached. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., September 9 , 1948. 


APPENDIX A 

THE EFFECT OF SWEEP ON THE FREQUENCIES OF A CANTILEVER BEAM 


Early in the investigation it was decided to make an experi- 
mental vibration study of a simple beam at various sweep 
angles. The uniform, plate-like aluminum-alloy beam shown 
in figure 1 8 was used to make the study amenable to analysis. 
Length-chord ratios of 6, 3, and 1.5 were tested, the length / 
being defined as the length along the midchord. A single 60- 
inch beam was used throughout the investigation, the desired 
lengt h and sweep angle being obtained by clamping the beam 
in the proper position with a 1 %- by 1 !4- by 14-incli 
aluminum-alloy crossbar. 

Figures 18 and 19 show the variation in modes and fre- 
quencies with sweep angle. In most cases, an increase in 
sweep angle increased the natural vibrational frequencies. 
As expected, the effect of sweep was more pronounced at the 
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Figure 18— Change in nodal lines with sweep and length-chord ratio for the vibration o 
an aluminum-alloy beam. 


smaller values of length-chord ratio. The fundamental mode 
was found by striking the beam and measuring the frequency 
with a self-generating vibration pick-up and paper recorder. 
The second and third modes were excited by light-weight 
electromagnetic shakers clamped to the beam. r I hese shakers 
were attached as close to the root as possible to give a node 
either predominantly spanwise or chordwise. The mode 
with the spanwise node, designated second mode, was pri- 
marily torsional vibration, whereas the mode with the chord- 
wise node, designated third mode, was primarily a second 
bending vibration. 

The first two bending frequencies and the lowest torsion 
frequency, determined analytically for a straight uniform 
unswept beam, are plotted in figure 19. Good agreement 
exists with the experimental results for the length-chord 
ratios of 6 and 3, but for a ratio of 1.5 (length equal to 12 
inches and chord equal to 8 inches) less favorable agreement 
exists. This discrepancy may be attributed to the fact that 
the beam at the short length-chord ratio of 1.5 resembled 
more a plate than a beam and did not meet the theoretical 
assumptions of a perfectly rigid base and of simple-beam 
stress distributions. The data are valid for use in comparing 
the experimental frequencies of the beam when swept with the 
frequencies at zero sweep, which was the purpose of the test. 


Experimental Analytical l/c 
o Q, 6 



Figure 19.— Variation of frequencies with sweep and length-chord ratio for the vibration 
of an aluminum-alloy beam. 
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APPENDIX B 


DISCUSSION OF THE REFERENCE FLUTTER SPEED 


For use in comparing data of swept and unswept wings, a 
reference flutter speed V R is convenient. This reference 
flutter speed is the flutter speed determined from the simpli- 
fied theory of reference 7. This theory deals with two- 
dimensional unswept wings in incompressible flow and de- 
pends upon a number of wing parameters. The calculations 
in this report utilize parameters of sections perpendicular to 
the leading edge, first bending frequency, uncoupled torsion 
frequency, density of testing medium at time of flutter, and 
zero damping. Symbolically, 


Since m is proportional to 6 , 


Ttpu . 

k= = k 0 cos A 

m 


Similarly, since / is proportional to b , 

, 0.56 lEI fr\, . , 

A = ~ 12 ~ (cos A) 2 

Also, because j a is independent of A, 


1 R b(x3 a f ^/C, X C g, X ea , V a" , f ^ 


A 

fa 



(cos A) 2 


Variation in reference flutter speed with sweep angle for 
sheared swept wings. — The reference flutter speed is 
independent of sweep angle for a homogeneous rotated wing 
and for homogeneous wings swept back by keeping the length- 
chord ratio constant. For a series of homogeneous wings 
swept back by the method of shearing, however, a definite 
variation in reference flutter speed with sweep angle exists 
since sweeping a wing by shearing causes a reduction in 
chord perpendicular to the wing leading edge and an increase 
in length along the mid chord as the angle of sweep is in- 
creased. The resulting reduction in the mass-density-ratio 
parameter and first bending frequency tends to raise the 
reference flutter speed, whereas the reduction in semichord 
tends to lower the reference flutter speed as the angle of 
sweep is increased. The final effect upon the reference 
flutter speed depends on the other properties of the wing. 
The purpose of this section is to show the effect of these 
changes on the magnitude of the reference flutter speed for 
a series of homogeneous sheared wings having properties 
similar to those of the sheared swept models used in this 
report. 

Let the subscript 0 refer to properties of the wing at zero 
sweep angle. The following parameters are then functions 
of the sweep angle: 


6 = 6 0 cos A 



cos 


A 


An estimate of the effect on the flutter speed of these 
changes in semichord and mass parameter with sweep angle 
may be obtained from the approximate formula given in 
reference 7, 


1 r~ be 


0.5 


0.5- 


a-\-x a ~ ^ /fi ° ' v COS A 


1 his approximate analysis of the effect on the reference 
flutter speed does not depend upon the first bending frequency 
but assumes f h /f a to be small. 

In order to include the effect of changes in bending- torsion 
frequency ratio, a more complete analysis must be carried 
out. Figure 20 presents some results of a numerical analysis 
based on a homogeneous wing with properties at zero sweep 
angle as follows: 


r c *=50 
X e a — 4 5 
r a 2 =0.2 5 
f a = 100 

In figure 20 the cur ve show ing the decrease in V R with A is 
slightly above the V cos A factor indicated by the approxi- 
mate formula. 

Effect of elastic-axis position on reference flutter speed. — 

As pointed out in the definition of elastic axis, the measured 


6o=0.3:i3 
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Sweep angle , A, deg 


Figure 20. — ' Variation in reference flutter speed with sweep for sheared wings. 


locus of elastic centers x ea ' fell behind the section elastic axis 
for the swept models with bases parallel to the air stream. 
In order to get an idea of the effect of elastic-axis position on 
the chosen reference flutter speed, computations were made 
both of V R and a second reference flutter speed V R ' similar 
to V R except that x ea ' was used in place of x ea . The maximum 
difference between these two values of reference flutter speed 
was of the order of 7 percent. This difference occurred at a 
sweep angle of 60° when the wing elastic axis was farthest 
behind the section elastic axis. Thus, for wings of this type, 


the reference flutter speed is not very sensitive to elastic-axis 
position. The reference flutter frequency j R was found in 
conjunction with V R . The maximum difference between 
j R and ji i was less than 10 percent. Thus, the convenient 
use of the reference flutter speed and reference frequency is 
not altered by these elastic-axis considerations. 
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TABLE I.— DATA FOR SHEARED 


Model 

A 

(deg) 

A t 

Aj 

(cps) 

Aj 

(cps) 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

M tt 

l 

(in.) 

C 

(in.) 

b 

(ft) 

X" 

(ix'reent 

chord) 

X t a 

(percent 

chord) 

lea' 

(percent 

chord) 

a+Ta 

a 

r a 2 

1 

K 

P 

/slugs \ 
\cu ft/ 

Percent 
Freon ; 

Spruce wings 

11 A 

0 

2 

45 

272 

108 

107 

15, 000 

25, 100 

16-005 

0. 89 

16.0 

8.0 

0. 333 

48.4 

45 

45 

-0. 032 

-0. 10 

0. 232 

13.3 

0. 00287 

95 

11 A' 

0 

2 

26 



59 

37 





16-005 

.89 

16.0 

8.0 

. 333 

48.4 

26.6 

26.6 

-. 032 

-. 468 

.396 

17.6 

.00217 

0 

11B' 

0 

2 

20 


61 

43 





16-005 

.89 

16.0 

8.0 

. 333 

48.4 

29.7 

29.7 

-.032 

-.406 

.371 

40.5 

. 000943 

88 

12 

15 

2 

43 



103 

103 

14,400 

54, 700 

16-005. 2 


16. 6 

7. 72 

.321 

48.5 

46.3 

46 

-. 03 

-.074 

.23 

5. 69 

. 00725 

96 

12 

15 

2 

42 


105 

105 

14,400 

54 700 

16-005. 2 

.88 

16. 6 

7. 72 

.321 

48.5 

46.3 

46 

-.03 

-.074 

. 23 

8.47 

. 00486 

98 

12 

15 

2 

42 


103 

102 

14,400 

54, 700 

16-005. 2 

.88 

16. 6 

7. 72 

.321 

48.5 

46.3 

46 

-.03 

-.074 

. 23 

11. 2 
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97 
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30 

2 

33 

196 

94 

93 

11, 100 

53, 500 

16-005.8 

.87 

IS. 2 

6. 87 

. 284 

48.8 

46.0 

49 

-. 024 

-.080 

.23 

7. 15 

. 00746 

99 

14 

45 

2 

22 

139 

93 

92 

9, 240 

33, 000 

16-007. 1 

.85 

22.6 

5.62 

. 23 1 

48.8 

46. 0 

60 

-. 024 

-.080 

23 

7. 78 

. 00720 

85 

14 

45 

2 

21 

136 

92 

91 

9, 240 

33, 000 

16-007. 1 

.85 

22. 6 

5. 62 

. 234 

48.8 

46. 0 

60 

-.024 

-.080 

.23 

19.8 

. 00285 

94 

15 

00 

2 

12 

68 

93 

93 

4, 520 

19, 100 

16-010 

.81 

32. 0 

4.0 

. 1 67 

48.8 

46.0 

65 

-. 024 

-.080 

.23 

9. 10 

. 00757 

92 

15 

60 

2 

12 

67 

93 

93 

4, 520 

19, 100 

16-010 

.81 

32. 0 

4.0 

. 167 

48.8 

46.0 

65 

-.024 

-.080 

.23 

14.0 

. 00493 

90 


Balsa wings 


22' 

15 

2 

31 

155 

63 

61 



16-005. 2 

0.88 

16.6 

7. 72 

0.321 

48.8 

42.4 

42.4 

-0. 024 

-0. 152 

0. 292 

2. 19 

0. 00854 

98 

22' 

15 

2 

31 

154 

64 

62 



16-005. 2 

.88 

16.6 

7. 72 

. 321 

48.8 

42.4 

42.4 

-.024 

-. 152 

.292 

3. 82 

. 00488 

93 

22' 

15 

2 

31 

154 

64 
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.88 

16. 6 

7. 72 

.321 

48.8 

42.4 

42.4 

-.024 

-. 152 

. 292 

18.7 

. 00100 

92 

23 

30 

2 

35 

219 

89 

89 

6, 230 

27, 900 

16-005. 8 

.87 

18.2 

6. 87 

.284 

48. 0 

48. 0 

52 

-.04 

-.04 

.304 

3. 18 

. 00864 

99 

23 

30 

2 

34 

216 

89 

89 

6. 230 

27. 900 

16-005.8 

.87 

18.2 

6. 87 

.284 

48. 0 

48. 0 

52 

-.04 

-.04 

. 304 

8. 54 

.00321 

91 

23 

30 

2 

34 

220 

91 

91 

6, 230 

27, 900 

16-005.8 

.87 

18.2 

6. 87 

. 284 

48.0 

48.0 

52 

-.04 

-.04 

.304 

9. 15 

. 00300 

89 

23 

30 

2 

34 

216 

89 

89 

6. 230 

27. 900 

16-005.8 

.87 

18. 2 

6. s7 

.284 

48.0 

48. 0 

52 

-.04 

-.04 

. 304 

14.9 

. 00184 

90 

24 

45 

2 

19 
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73 

73 

2,810 

10. HID 

16 007. 1 

.85 

21.8 

5. 66 

.236 

47.0 

49.0 

57 

-.06 

-.02 

.311 

3. 64 

. 00784 

85 

24 

45 

2 

19 
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75 

75 

2,810 

10, 800 

16-007. 1 

.85 

21. s 

5.66 

. 236 

47.0 

49.0 

57 

-.06 

-.02 

.311 

8 . 40 

. 00339 

93 

24 

45 

2 

19 

122 

75 

75 

2,810 

10. 800 

16-007. 1 

.85 

21.8 

5. 66 

. 236 

47.0 

49.0 

57 

-.06 

-.02 

.311 

13.2 

.00216 

91 

24 

45 

2 

19 

120 

74 

74 

2,810 

10, 800 

16-007. 1 

.85 

21.8 

5. 66 

. 236 

47. 0 

49.0 

57 

-. 06 

-.02 

.311 

29.4 

. 000970 

74 

24 

45 

2 

19 

120 

73 

73 

2,810 

10. 800 

16-007. 1 

.85 

21.8 

5. 66 

.236 

47.0 

49.0 

57 

-.06 

-.02 

.311 

30.6 

.000933 

89 

25 A 

60 

2 

8.6 

54 

66 

65 

1,950 

6, 470 

16-010 

.81 

32.0 

4.0 

. 167 

46.9 

40. 0 

71 

-. 062 

-. 20 

. 359 

34.6 

. 000954 

88 

25 B 

60 

2 

8.6 

48 

70 

68 


5,500 

16-010 

.81 

32.0 

4.0 

. 167 

46.9 

40.0 

71 

-.062 

-.20 

.359 

9.36 

. 00353 

91 


TABLE II— ROTATED 


Model 

A 

(deg) 

A, 

A, 

(cps) 

Aj 

(cps) 

ft 

(cps) 

A 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

M er 

/ 

(in.) 

c 

(in.) 

b 

(ft) 

X eg 

(percent 

chord) 

Xta 

(percent 

chord) 

(percent 

chord) 

d+Xa 

a 

Ta 2 

K 

P 

/slugs\ 
\cu ft/ 

Percent 

Freon 

Lengthwise laminations 

30A 

0 

6.20 

11. 9 

76. 0 

90.4 

83. 0 

3, 760^. 


16-010 

0.81 

24.8 

4 

0. 167 

46.0 

35 

35 

-0.08 

-0. 30 

0.311 

36.8 

0. 00220 

0 

30 B 

0 

6. 20 

12.0 

72.6 

90.0 

88.0 

3, 760 

6, 920 

16-010 

.81 

24.8 

4 

. 167 

46.0 

40 

40 

-.08 

-.20 

.277 

37.8 

.00214 

0 

30 B 

30 

4.65 

12. 1 

73. 0 

91.0 

88.8 

3, 760 

6, 920 

16-010 

.81 

24.8 

4 

. 167 

46.0 

40 

40 

-.08 

-.20 

.277 

37. 7 

. 00215 

0 

30 B 

30 

4.65 

12.0 

73. 0 

90. 0 

88.0 

3, 760 

6. 920 

16 <>iii 

.81 

24.8 

4 

. 167 

46. 0 

40 

40 

-.08 

-.20 

.277 

37.8 

.00214 

0 

30 B 

45 

3. 10 

12. 1 

73. II 

91.0 

88.8 

760 

6. 920 

16-010 

.81 

24.8 

4 

. 167 

46.0 

40 

40 

-.08 

-. 20 

. 277 

37. 8 

. 00214 

0 

30 B 

45 

3. 10 

12.2 

73.0 

90. 0 

88.0 

3, 760 

6, 920 

16-010 

.81 

24.8 

4 

. 167 

46.0 

40 

40 

-.08 

-. 20 

.277 

37.8 

.00214 

0 

30 B 

60 

1.55 

12.0 

72.5 

90.0 

88.0 

3, 760 

6,920 

16 Dili 

.81 

24.8 

4 

. 167 

46.0 

40 

40 

-.08 

-. 20 

.277 

39.8 

. 00204 

0 

30C 

0 

6. 20 

12.2 

69. 0 

86. 0 

75.8 

4. 000 

6. 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-. 22 

.292 

40.5 

. 00200 

89 

30C 

0 

6. 20 

12.2 

69. 0 

86. 0 

75.8 

4, (KM) 

6, 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

.292 

98.9 

. 0(X)S20 

86 

30C 

0 

6. 20 

13.3 

70. 0 

84.0 

74.2 

4, 000 

6. 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

. 292 

92.6 

. 000876 

83 

30C 

15 

5. 78 

12.2 

69. 0 

86. 0 

75.8 

4, 000 

6, 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

.292 

92.6 

. 000870 

81 

30C 

30 

4.65 

12.2 

69.0 

86.0 

75. 8 

4, 000 

6, 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

.292 

40.0 

. 00202 

89 

30C 

30 

4.65 

12.2 

70.0 

86.5 

76.2 

4,000 

9, 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

.292 

81.4 

. 000995 

86 

30C 

30 

4. 65 

12.2 

70. 0 

86. 5 

76.2 

4, 000 

6, 950 

16-010 

.81 

24.8 

4 

. 167 

48. 5 

39 

39 

-. 03 

-.22 

.292 

80.0 

. 00100 

85 

30C 

45 

3. 10 

12.2 

70. 0 

86.5 

76.2 

4. (M)0 

6. 950 

16-010 

.81 

24.8 

4 

. 167 

48.5 

39 

39 

-.03 

-.22 

.292 

45.2 

.00179 

87 

301) 

15 

5. 78 

13.2 

80.2 

87. 1 

82.4 

4, 350 


16 OK) 

.81 

24.8 

4 

. 167 

48 

39.5 

39. 5 

-.04 

-.21 

.280 

8. 70 

. 00933 

99 

30D 

15 

5. 78 

13.2 

80.2 

87. 1 

82.4 

4, 350 


16-010 

.81 

24.8 

4 

. 167 

48 

39. 5 

39.5 

-.04 

-.21 

.280 

8. 72 

. 00930 

99 

301) 

15 

5.78 

13.2 

80.2 

87. 1 

82.4 

4, 350 


16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39.5 

-.04 

-.21 

. 280 

8. 76 

. 00927 

99 

301) 

30 

4.65 

13.5 

81. 7 

92. 5 

87.4 

4, 350 


16 OKI 

.81 

24.8 

4 

. 167 

48 

39. 5 

39.5 

-.04 

-.21 

. 280 

8. 90 

.00910 

99 

301) 

45 

3. 10 

13.3 

81. 7 

88. 2 

83. 4 

4, 350 


16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39. 5 

-.04 

-.21 

. 280 

8. 85 

. 00905 

99 

30 D 

60 

1.55 

13.5 

82.0 

90.5 

85.5 

4, 350 


16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39. 5 

-.04 

-.21 

.280 

9. 54 

. 00852 

99 

Chord wise laminations 

40 A 

0 

6. 20 

9.4 

57.4 

90.0 

88.4 

3, 540 

5, 250 

16-010 

0.81 

24.8 

4 

0. 167 

46 

40 

40 

-0. 08 

-0. 20 

0. 277 

36.5 

0. 00222 

0 

40A 

0 

6.20 

9. 6 

57. 1 

91.0 

88.5 

3, 540 

5. 250 

16-010 

.81 

24.8 

4 

. 167 

46 

40 

40 

-.08 

-.20 

.277 

24.2 

. 00334 

90 

40 A 

0 

6. 20 

9.6 

57. 1 

91.0 

88.5 

3, 540 

5, 250 

16-010 

.81 

24.8 

4 

. 167 

46 

40 

40 

-.08 

-.20 

.277 

37. 7 

.00215 

89 

40 A 

0 

6. 20 

9.6 

57. 1 

91.0 

88.5 

3, 540 

5, 250 

16-010 

.81 

24.8 

4 

. 167 

46 

40 

40 

-.08 

-.20 

.277 

75. 0 

. 00108 

82 

40 A 

15 

5. 78 

9.3 

55.8 

90. 6 

88.2 

3, 540 

5, 250 

16-010 

.81 

24.8 

4 

. 167 

46 

40 

40 

-.08 

-.20 

.277 

35. 1 

.00231 

0 

40A 

30 

4.65 

9. 3 

55. 8 

90.6 

88.2 

3, 540 

5. 250 

16 010 

.81 

24.8 

4 

. 167 

46 

40 

40 

-.08 

-.20 

.277 

37.5 

.00216 

0 

40B 

0 

6. 20 

9.5 

55.0 

90. 5 

85.5 

3,710 

5, 020 

16-010 

.81 

24.8 

4 

. 167 

49 

40 

40 

-.02 

-.20 

.297 

35.5 

.00228 

0 

40C 

0 

6. 20 

9.0 

54.4 

61.0 

58. 2 

2, 280 

4,350 

16-010 

.81 

24.8 

4 

. 167 

46 

38.5 

38. 5 

-.08 

-.23 

.287 

8. 74 

. 00928 

100 

40D 

0 

6. 20 

9.4 

58.0 

88.9 

84.0 

3. 330 

5. 050 

16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39. 5 

-.04 

-.21 

.280 

79. 0 

. 000969 

84 

401) 

15 

5.78 

9.6 

58.3 

88.9 

84.0 

3, 330 

5, 050 

16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39. 5 

-.04 

-.21 

.280 

36.2 

.00212 

89 

401) 

15 

5. 78 

9.5 

57.9 

87.5 

82.6 

3, 330 

5, 050 

16-010 

.81 

24.8 

4 

. 167 

48 

39. 5 

39.5 

-.04 

-.21 

.280 

80.0 

. 000956 

87 

401) 

30 

4.65 

9.5 

57.5 

89.0 

84. 1 

3. 330 

5. 050 

16-010 

.81 

24.8 

4 

. 167 

48 

39.5 

39.5 

-.04 

-.21 

.280 

88.2 

.0 him;: 

85 

40 D 

45 

3. 10 

9.6 

58.3 

88.9 

84.0 

3, 330 

5. 050 

16-010 

.81 

24.8 

4 

. 167 

48 

39. 5 

39.5 

-.04 

-.21 

.280 

39. 1 

. 00196 

86 


TABLE III— DATA FOR MODELS USED IN 


Model 

A 

(deg) 

A t 

A, 

(cps) 

Aj 

(cps) 

(cps) 

A 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

M er 

l 

(in.) 

c 

(in.) 

b 

(ft) 

Xeg 

(percent 

chord) 

Xea 

(percent 

chord) 

X ea ' 

(percent 

chord) 

a+.r a 

a 

Ta~ 

L 

P 

/ slugs \ 
Veil ft / 

Percent 

Freon 

50 A 

-30 

4.65 

15 

87 

168 

137 

10, 100 

14, 100 

16-010 

0.81 

24.8 

4 

0. 167 

50 

33 

33 

0.0 

-0.34 

0. 352 

7. 98 

0. 00895 

96 

50 A 

-15 

5. 78 

15 

87 

168 

137 

10. 100 

14, 100 

16-010 

.81 

24.8 

4 

. 167 

50 

33 

33 

.0 

-.34 

.352 

8. 00 

.00892 

96 

50 A 

0 

6.20 

15 

87 

ia3 

133 

10, 100 

14, 100 

16-010 

.81 

24.8 

4 

.167 

50 

33 

33 

.0 

-.34 

.352 

33.1 

.00216 

0 

50B 

0 

6.20 

14 

82 

166 

116 

11,400 

11,900 

16 010 

.81 

24.8 

4 

. 167 

50 

26 

26 

.0 

-.48 

. 456 

8.66 

.00823 

99 

50B 

15 

5. 78 

14 

80 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

.167 

50 

26 

26 

.0 

-.48 

. 456 

8. 58 

.00831 

99 

50 B 

30 

4.65 

14 

80 

166 

116 

11,400 

11,900 

16-010 

.81 

24.8 

4 

. 167 

.50 

26 

26 

.0 

-.48 

.456 

9. 04 

.00787 

99 

50 B 

45 

3. 10 

14 

80 

166 

116 

1 1 , 400 

11,900 

16-010 

.81 

24.8 

4 

.167 

.50 

26 

26 

.0 

-.48 

.456 

9.45 

. 00756 

99 
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SWEPT MODELS— SERIES I 


(/ps) 

fa 

(cps) 

/a 

(cps) 

f. 

fa 

U 

/r 

f. 

U 

<P 

(deg) 

(si) 

M 

V e 

(mph) 

Vr 

(mph) 

Vr’ 

(mph) 

Fa 

(mph) 

v , 

blt)a 

Vs 

Vr 

Fa 

Vd 

(mph) 

Remarks 

Spruce wings 

06 

42 

38 

64 

62 

55 

61 

54 

37 

37 

36 

70 

40 

42 

70 

71 
69 
60 
56 
51 
53 
51 

71 

69 

65 

61 

53 

58 

55 

0. 62 
1 . 12 
.87 
.63 
.59 
.54 
.66 
.59 
.41 
.40 
.39 

0. 93 
1.03 
.91 
.92 
.87 
.80 
1.01 
.97 
.72 
.70 
.71 

0. 87 
.80 
.94 
.88 
.70 
.64 
.65 

50 

80 

170 

70 

50 

50 

70 

60 

40 

40 

30 

235 

85.0 

70.5 

375 

320 

307 

334 

300 

234 

265 

264 

0. 82 
.24 
.74 
.64 
.71 
.79 
.62 
.56 
.81 
.51 
.62 

274 
191 
262 
218 
245 
276 
202 
196 

275 
179 
222 

260 

129 

197 

176 

206 

225 

154 

134 

191 

103 

124 

260 

129 

197 

107 ’ 
127 

205 

220 

159 

166 

245 

184 

222 

1.80 
3. 58 
4.22 
1.54 
1.70 
1.95 
1.77 
2. 11 
2. 99 
2. 71 
3.35 

1.05 
1.48 
1.33 
1.24 
1. 19 
1. 23 
1.31 
1.46 
1.44 
1.73 
1.79 

1.20 
1.25 
1.27 
1.18 
1. 12 
.97 
1.00 

314 

583 

183 

175 

217 

245 

149 

119 

187 

105 

122 

Tunnel excitation frequency, 67 cps. 

YT,„D>1 fciiarf / Slotted AH inches from trailing 
Model failed .| edge slots uncovered. 

Tunnel excitation frequency, 61 cps. 

Balsa wings 

50 

51 
45 
60 
62 
60 
53 
51 
49 
45 

34" 

29 

46 

48 
46 
62 
62 
63 
60 

49 
49 
48 

44 
43 
37 

45 

48 

46 

65" 

65 

58 

’45" 

40 

48 

0. 82 
.83 
.72 
.68 
.70 
.67 
.59 
. 71 
.65 
.60 

.47 

.44 

1.07 

1.07 

.96 

.96 

1.01 

.96 

.87 

1.06 

1.00 

.95 

".’79 

.75 

1. 06 
.98 

’.’95 

‘."82 

’.’84 

. 75 
.72 

70 

50 

50 

130 

70 

60 

90 

90 

40 

40 

60 

10 

101 

74.7 
54.2 

189 

152 

171 

152 

125 

120 

108 

83.5 
78.0 

76.8 

73.6 

0. 30 
.34 
.64 
.42 
.62 
.66 
.81 
.34 
.54 
.64 
.76 
.81 
.79 
.41 

1 

104 

119 

224 

142 

212 

229 

275 

121 

180 

215 

281 

277 

272 

139 

97.3 

95.0 
167 
137 
176 
185 
221 

97.1 
132 
160 
226 
226 
161 

93.5 

169 

97.5 

96 

168 

176 

22 1 

145 

252 

278 

161 

1.25 
1.41 
2. 64 
1.31 
1.95 

2. 07 
2.53 
1.63 
2.35 
2. 82 

3. 76 
3. 77 
5. 90 
2. 85 

1.07 

1.25 

1.34 
1.04 
1.21 
1.24 

1.24 

1.25 
1.37 

1.35 
1 . 25 
1 . 22 
1.69 
1.49 

1. 24 
1.33 

1.21 

1.23 

I. 24 

i.'io 

.98 

.86 

79.9 
107 
238 
110 
180 
190 
237 
80. 1 
127 
159 
232 
232 
210 
115 

1 Tunnel excitation frequency, 49 cps. 
j Slotted 2Mfi inches from trailing edge. 

Tunnel excitation frequency, 61 cps. 
Tunnel excitation frequency, 61 cps. 

Model failed. 

Model failed. 


WINGS— SERIES II 


fe 

(cps) 

/« 

(cps) 

/a 

(cps) 

fe 

fa 

fe 

fH 

fe 

/A 

<P 

(deg) 

(i) 

M 

V e 

(mph) 

Vr 

(mph) 

Vr' 

(mph) 

Fa 

(mph) 

v e 

bo)a 

Ve 

Vr 

Ve 

Fa 

V d 

(mph) 

Remarks 

Lengthwise laminations 

42 

45 


0.51 

0.91 


70 

127 

0. 30 

232 

209 

209 


3.91 

1. 11 


318 

Wing failed; tunnel excitation frequency, 















40.7 cps. 

48 

44 

46 

.54 

1.08 

1.04 

60 

121 

.29 

229 

212 

212 

215 

3. 64 

1.08 

1.06 

263 


51 

47 

47 

.57 

1.08 

1.08 

60 

126 

.30 

235 

214 

214 

229 

3. 74 

1. Ki 

1.03 

266 


50 

44 

47 

. 57 

1. 14 

1.08 

40 

129 

.30 

237 

212 

212 

229 

3. 77 

1. 12 

1.04 

263 


44 

47 





166 

.34 

269 

214 

214 

265 

4. 28 

1.26 

1.01 

266 


’55 

44 

47 

.62 

1. 25 

Lie 

’50" 

169 

.35 

272 

212 

212 

265 

4. 32 

1.28 

1.02 

263 

Wing failed. 

46 

48 




275 

.45 

350 

219 

219 

353 

5. 59 

1.60 

.99 

265 

34 

41 

.*45 

.’83 


’36’ 

104 

.63 

219 

189 

189 


4. 05 

1. 16 



249 


24 

37 


.32 

. 66 


30 

74.4 

.81 

286 

290 

290 


5. 29 

.986 



393 

Wing failed. 

21 

36 


.29 

.59 


30 

79. 6 

.82 

288 

270 

270 


5. 43 

1.07 


369 

27 

36 


.36 

.74 


30 

72. 5 

.78 

278 

282 

282 


5. 13 

.986 



376 


37 

41 


.48 

.89 


50 

113 

.65 

226 

187 

187 


4. 18 

1.21 


248 


41 






88. 1 

.81 

284 

263 

263 


5. 22 

1.08 



355 


31 

38 


.'40 

.’80 


§5 

88.6 

.81 

289 

260 

260 


5. 32 

1. 11 



352 


40 

41 


.53 

.98 


40 

147 

.76 

273 

199 

199 


5. 02 

1.37 

1. 04 

265 


50 

51 

‘51 

.61 

.98 

.98 

50 

110 

.31 

104 

100 

100 

166 

1. 77 

1.05 

119 


51 

52 

52 

.61 

.98 

. 98 

50 

115 

.32 

107 

100 

100 

101 

1.82 

1.08 

1.06 

119 


51 

52 

52 

.61 

.98 

.98 

50 

121 

.33 

109 

1(10 

100 

101 

1.85 

1. 10 

1.08 

119 


53 

54 

56 

.61 

.98 

.95 

40 

150 

.38 

123 

106 

106 

116 

1.97 

1. 16 

1.06 

129 


56 

52 

55 

.67 

1.08 

1.02 

60 

178 

.41 

135 

101 

101 

130 

2. 26 

1.34 

1.04 

122 


65 

53 

57 

.77 

1.24 

1. 14 

90 

307 

.55 

182 

107 

107 

182 

2. 98 

1.70 

1.00 

130 


Chord wise laminations 

62 

47 


0. 70 

1.33 


140 

82. 0 

0. 24 

188 

211 

211 


2. 98 

0. 892 


260 

Tunnel excitation frequency, 57 cps. 

56 

49 


.63 

1. 15 


60 

86. 7 

. 45 

155 

184 

184 


2. 45 

.843 



212 

61 

46 


. 69 

1.33 


70 

69.2 

.50 

172 

215 

■J 1 5 


2. 72 

.800 


265 


61 

43 


. 69 

1.44 


70 

63. 6 

.65 

234 

299 

299 


3. 70 

. 784 



373 


61 

46 


.68 

1.30 


90 

93. 9 

.26 

201 

208 

208 


3. 19 

.967 



254 

Wing failed. 

46 






127 

.30 

235 

213 

213 


3. 73 

1.10 



263 

61 

45 


.71 

1.37 


"io" 

77.7 

.23 

178 

191 

191 


2.91 

.932 



247 

Wing failed. 

29 

36 


. 51 

.83 


80 

57.6 

.23 

75.3 

74. 5 

74.5 


1.81 

1.01 



90. 1 

Wing failed. 

62 

40 


.73 

1.54 


30 

52.3 

.62 

221 

281 

281 


3. 69 

.787 



370 

Tunnel excitation frequency, 61 cps. 

62 

44 


.74 

1.41 


70 

72.7 

.51 

177 

194 

194 


2. 95 

.913 



251 


61 

40 


.71 

1.54 


50 

57.9 

.67 

236 

279 

279 


3.99 

. 846 



367 


65 

40 


. 77 

1.63 


60 

79. 1 

.82 

290 

• 20 s 

298 


4.83 

.973 



392 

Wing failed. 

32 

44 

::: 

.38 

.73 

-- 

80 

138 

.73 

254 

200 

200 

.... 

4.24 

1.27 

— 

261 

iWEEPFORWARD TESTS— SERIES 

III 











(cps) 

/« 

(Cps) 

/A 

(cps) 

fe_ 

fa 

u 

fR 

fe 

/A 

<P 

(deg) 

(si) 

M 

V e 

(mph) 

Vr 

(mph) 

Vr' 

(mph) 

Fa 

(mph) 

Ve 

boi a 

Ve 

Vr 

Vs 

Fa 

V D 

(mph) 

Remarks 


98 

98 






73.4 

0. 26 

86.9 

174 

174 


0. 888 

0. 498 


294 

Model diverged. 







107 

.31 

105 

174 

174 


1.075 

. 603 


294 

Model diverged. 

102 

91 

84 

74 

79 


0. 77 
.78 
.72 
. 63 

1.29 

.97 

.90 

.80 


40 

100 

70 

180 

211 

260 

257 

352 

. 40 

303 

319 

319 


3. 18 

.949 


579 


94 

94 

93 

93 



.52 

.51 

.61 

170 

L69 

202 

172 

172 

172 

172 

179 


2. 05 

. 989 


704 





2. 04 

.982 


700 

Model failed. 



179 


2. 44 

1.125 


720 


98 


.84 

1.05 


100 

423 

.68 

226 

179 

179 


2.73 

1.265 


736 
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TABLE IV.— SWEPT MODELS OF A CONSTANT 


Model 

A 

(deg) 

A K 

(cps) 

/a 2 

(cps) 

ft 

(cps) 

/« 

(cps) 

GJ 

(lb-in. 2) 

El 

(lb-in. 2) 

NACA 

airfoil 

section 

Id" 

l 

(in.) 

(in.) 

b 

(ft) 

Xeg 

(percent 

chord) 

(percent 

chord) 

X ea 

(percent 

chord) 

a Xa 

a 

r a 2 

62 

15 

7.95 

4.9 

29.1 

72.5 

71.8 

3, 730 

7, 820 

16-010 

0.81 

34 

4 

0. 167 

41 

44 

46 

-0. 18 

-0. 12 

0. 175 

62 

15 

7. 95 

4.9 

29. 1 

73.4 

72.5 

3. 730 

7, 820 

16-010 

.81 

34 

4 

. 167 

41 

44 

46 

— . is 

— . 12 

. 175 

62 

15 

7. 95 

4.9 

29.1 

73.4 

72.5 

3. 730 

7, 820 

16-010 

.81 

34 

4 

. 167 

41 

44 

46 

-. 18 

-. 12 

. 175 

62 

15 

7. 9c 

4.9 

29. 6 

73.5 

72.7 

3,730 

7,820 

16-010 

.81 

34 

4 

. 167 

41 

44 

46 

— . 18 

— . 12 

. 175 

63 

30 

6. 38 

4.6 

25.8 

73. 5 

73.0 

5. 450 

5, 870 

16 OK) 

.81 

34 

4 

. 167 

41 

44 

47 

— . 18 

-. 12 

. 175 

63 

30 

6. 38 

3.9 

24.0 

73. 0 

72.4 

5. 450 

5, 870 

16-010 

.81 

34 

4 

. 167 

41 

44 

47 

-. 18 

-. 12 

. 175 

as 

30 

6. 38 

4.6 

25.8 

73. 5 

73.0 

5, 450 

5, 870 

16-010 

.81 

34 

4 

.167 

41 

44 

47 

-. 18 

-. 12 

. 175 

as 

30 

6. 38 

4. 0 

24.0 

73.0 

72.4 

5, 450 

5, 870 

16-010 

.81 

34 

4 

. 167 

41 

44 

47 

-. 18 

-. 12 

. 175 

as 

30 

6. 38 

4.0 

24.0 

73.0 

72.4 

5, 450 

5, 870 

16-010 

.81 

34 

4 

. 167 

41 

44 

47 

-. 18 

-. 12 

. 175 

64 

45 

4.75 

4.4 

29. 0 

66.0 

65.5 

3. 500 

6, 080 

16-010 

.81 

34 

4 

. 167 

41 

44 

57 

-. 18 

-. 12 

. 175 

64 

45 

4. 75 

4.2 

27.0 

66. 0 

65. 5 

3, 500 

6, 080 

16-010 

.81 

34 

4 

. 167 

41 

44 

57 

-. 18 

-. 12 

. 175 

64 

45 

4.75 

4.2 

27.0 

66. 0 

65.5 

3, 500 

6, 080 

16-010 

.81 

34 

4 

.167 

41 

44 

57 

— . 18 

-. 12 

. 175 

64 

45 

4.75 

4. 1 

27. 0 

65.0 

64.4 

3, 500 

6, 080 

16 010 

.81 

34 

4 

. 167 

41 

44 

57 

— . 18 

-. 12 

. 175 

64 

45 

4.75 

4. 1 

27.0 

65.0 

64.4 

3, 500 

6, 080 

16-010 

.81 

34 

4 

. 167 

41 

44 

57 

-. 18 

— . 12 

. 175 

6.5 

60 

2. 12 

5.7 

33. 4 

77.0 

76.2 

4. 650 

11.980 

16-010 

.81 

34 

4 

. 167 

41 

44 

71 

-. 18 

-. 12 

. 175 


TABLE V.— DATA FOR SWEPT MODELS OF A 


Model 

A 

(deg) 

A t 

/*, 

(cps) 

fK 

(cps) 

ft 

(cps) 

fa 

(cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

M r 

1 

(in.) 

C 

(in.) 

b 

(ft) 

X eg 

(percent 

chord) 

(percent 

chord) 

Xea 

(percent 

chord) 

a + .r a 

a 

T a“ 

72 

15 

6. 09 

7.6 

54 

97.3 

96.3 

3. 730 

7. 820 

16-010 

0.81 

26 

4 

0. 167 

41 

44 

46 

-0. 18 

-0. 12 

0. 175 

72 

15 

6. 09 

7. 6 

54 

97.3 

96.3 

3, 730 

7, 820 

16-610 

.81 

26 

4 

. 167 

41 

44 

46 

-. 18 

-. 12 

. 175 

73 

30 

4.88 

6.4 

40.0 

98.0 

97.0 

5, 450 

5, 870 

16-010 

81 

26 

4 

. 167 

41 

44 

47 

-. 18 

— . 12 

. 175 

73 

30 

4.88 

6.4 

40.0 

98. 0 

97.0 

5, 450 

5, 870 

16-010 

.81 

26 

4 

.167 

41 

44 

47 

-. 18 

-. 12 

. 175 

73 

30 

4.88 

6.4 

40.0 

98.0 

97.0 

5, 450 

5, 870 

16-610 

.81 

26 

4 

. 167 

41 

44 

47 

-. 18 

12 

. 175 

74 

45 

3. 25 

6.5 

40.0 

79.0 

78. 2 

3, 500 

6, 080 

16-010 

.81 

26 

4 

. 167 

41 

44 

57 

-. 18 

— . 12 

. 175 

74 

45 

3. 25 

6.7 

39. 5 

78.5 

77. 7 

3, 500 

6, 080 

16-610 

.81 

26 

4 

.167 

41 

44 

57 

-. 18 

— . 12 

. 175 

74 

45 

3.25 

6.7 

39. 5 

78.5 

77. 7 

3, 500 

6, 080 

16-610 

.81 

26 

4 

. 167 

41 

44 

57 

-. 18 

— . 12 

. 175 

75 

60 

1.65 

7.2 

51.8 

82.4 

81.6 

4. 650 

11.980 

16-010 

.81 

26 

4 

. 167 

41 

44 

71 

-. 18 

— . 12 

' 175 

75 

60 

1.65 

7.2 

51.8 

84.6 

as. 8 

4. 650 

11.980 

16-010 

.81 

26 

4 

. 167 

41 

44 

71 

-. 18 

-. 12 

! 175 


TABLE VI— DATA FOR TIP- 


Model 

A 

(deg) 

-1< 

/ a . 

(cps) 

/*, 

(cps) 

ft 

(cps) 

fa 

(cps) 

G.J 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

Me, 

/ 

(in.) 

(in.) 

b 

(ft) 

X eg 

(percent 

chord) 

X ea 

(percent 

chord) 

2 f tl 

(percent 

chord) 

a -1- Xa 

a 

r a 

84-1 

45 

3.fvl 

10 

60 

13S 

104 




16-610 

0.81 

29 

4 

0. 167 

51 

32 

44 

0. 02 

-0.36 

0.378 

84-2 

45 

3.63 

10 

61 

135 

107 



16-010 

.81 

29 

4 

.167 

51 

32 

44 

.02 

-.36 

.378 

84-3 

45 

3. as 

9.6 

58 

118 

93 



16-010 

.81 

29 

4 

.167 

51.5 

32 

44 

.03 

-.36 

.378 

85-1 

60 

2. 75 

5.0 

32 

92 

72 

10, 800 

13. 400 

16-010 

.81 

44 

4 

. 167 

50 

32 

58 

.0 

-.36 

.378 

85-2 

60 

2. 75 

5.0 

31 

95 

75 

9. 850 

12, 400 

16-010 

.81 

44 

4 

. 167 

50 

32 

58 

.0 

-.36 

.378 

85-3 

60 

2. 75 

5.0 

30 

80 

as 

11,200 

16.600 

16-010 

.81 

44 

4 

. 167 

51 

32 

58 

.02 

-.36 

.378 


TABLE VII. — DATA FOR MODELS USED TO DETERMINE 


Model 

A 

(deg) 

U 

(cps) 

/Aj 

(cps) 

ft 

(cps) 

fa 

(Cps) 

GJ 

(lb-in. 2 ) 

El 

(lb-in. 2 ) 

NACA 

airfoil 

section 

l 

(in.) 

(in.) 

b 

(ft) 

Xrg 

(percent 

chord) 

X ea 

(percent 

chord) 

i * -1 

(percent 

chord) 

a + Xa 

a 

Tor 

1 

91-1 

0 

6 

4.2 

24 

31 

23 

34, 100 

128, 000 

16-010 

48 

8 

0. 333 

29.9 

48 

48 

-0. 402 

-0. 04 

0. 307 

17.3 

91-2 

0 

6 

5.5 

36 

43 

43 

41.200 

108. 300 

16-010 

48 

8 

. 333 

41.0 

43.8 

43.8 

-.18 

-. 124 

. 179 

41 7 

91-2 

0 

6 

5.5 

36 

43 

43 

41.200 

ins. 300 

16-010 

48 

8 

. 333 

41.0 

43.8 

43.8 

-. 18 

-. 124 

. 179 

56. 4 

91-2 

0 

6 

5.3 

33 

42 

42 

41.200 

108. 3(H) 

16-010 

48 

8 

. 333 

41.0 

43.8 

43.8 

-. 18 

-. 124 

. 179 

12. 8 

91-2 

0 

6 

5.5 

36 

43 

43 

41.200 

108. 300 

16-010 

48 

8 

. 333 

41.0 

43.8 

43.8 

-. 18 

-. 121 

. 179 

95. 5 

91-3 

0 

6 

5.0 

30 

40 

40 

28. 500 

ai. 700 

16-010 

48 

8 

.333 

49. 0 

48.4 

48.4 

-.02 

-.032 

. 160 

44. 3 

91-3 

0 

6 

4.7 

29 

39 

39 

28. ,500 

83, 7i in 

16-010 

48 

8 

. 333 

49. 0 

48.4 

48.4 

-.02 

-.032 

. 160 

36. 4 

91-3 

0 

6 

4.7 

29 

39 

39 

28. 500 

83. 700 

16-010 

48 

8 

. 333 

49. 0 

48.4 

48.4 

-.02 

-.032 

. 160 

48. 4 

92-1 

15 

6. 09 

8.3 

48 

70 

62 

3, 730 

7. 820 

Modified 

16-010 

26 

4 

. 167 

31.2 

44 

46 

-.376 

-. 12 

.298 

77.9 

92-2 

15 

6. 09 

8.3 

49 

95 

95 

3, 730 

7, 820 

Modified 
16 OKI 

26 

4 

. 167 

42. 9 

44 

46 

-. 142 

-. 12 

. 136 

76.0 

92-3 

15 

6. 09 

8. 1 

47 

55 

52 

3, 730 

7, 820 

Modified 

16-010 

26 

4 

. 167 

54.5 

44 

46 

.090 

-.12 

.411 

74.5 

93-1 

30 

4.42 

6.3 

40 

78 

68 

5, 450 

5, 870 

Modified 

16-010 

23.6 

4 

. 167 

30 

44 

47 

-.40 

-. 12 

.310 

78.0 

93-2 

30 

4.42 

6.8 

44 

99 

99 

5, 450 

5, 870 

Modified 

16-010 

23. 6 

4 

. 167 

43 

44 

47 

-. 16 

-. 12 

.134 

74.0 

93-3 

30 

4.42 

6.3 

51 

54 

.50 

5. 450 

5. 870 

Modified 

16-010 

23. 6 

4 

. 167 

56 

44 

47 

. 12 

-. 12 

.428 

73.2 

94-1 

-(-45) 

3.81 

4.5 

26 

38 

35 

2. 120 

4. 520 

Modified 

16-010 

30. 5 

4 

. 167 

44.5 

56 


-. 11 

. 12 

.427 

68.2 

94-2 

-(-45) 

3.81 

4.8 

28 

70 

70 

2. 120 

4. 520 

Modified 
16 oio 

30.5 

4 

. 167 

57.0 

56 


. 14 

. 12 

. 134 

68.2 

94-3 

-(-45) 

3.81 

4.6 

28 

40 

38 

2. 120 

4. 520 

Modified 

16-010 

30.5 

4 

. 167 

69. 3 

.56 


.386 

. 12 

.307 

68.2 

95'- 1 

60 

1.65 

5. 6 


54 

50 

1.9(H) 

4.560 

Modified 
16 010 

26.4 

4 

.167 

31.4 

22 

41 

-.372 

-.56 

.267 

75.8 

95'-2 

60 

1.65 

5. 9 


71 

47 

1.900 

4. 560 

Modified 

16-010 

26.4 

4 

.167 

42.8 

22 

41 

-. 144 

-.56 

.308 

73.0 

95'-3 

60 

1.65 

5.8 

35 

40 

27 

1.900 

4, 560 

Modified 

16-010 

26.4 

4 

. 167 

54. 3 

22 

41 

.086 

-.56 

.779 

69. 0 
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LENGTH-CHORD RATIO OF 8.5— SERIES IV 


1 

P 

/ slugs \ 
V cu ft / 

Percent 

Freon 

(cps) 

/ft 

(cps) 

/. 

/« 

A 

/ft 

(deg) 

<7 

\ sq ft ) 

M 

Ve 

(mph) 

TV 

(mph) 

TV 

(mph) 

Ve 

tjUa 

T r , 
1 V 

tv 

(mph) 

Remarks 

13.5 

0. 00925 

99 

22 

35 

0.28 

0. 59 

30 

91.8 

0.29 

95.4 

105 

104 

1.85 

0. 905 

91.6 


37.6 

. 00333 

88 

20 

32 

.28 

.64 

20 

73.7 

.41 

143 

167 

171 

2. 76 

.856 

153 


59.5 

. 00210 

87 

19 

31 

.26 

.60 

20 

69.7 

.49 

175 

206 


3. 37 

. 850 

192 


130. 0 

.000964 

85 

16 

29 

.22 

.55 

20 

57.5 

.66 

234 

300 


4.50 

.780 

284 


15.2 

.00745 

73 

19 

35 

.27 

.56 

180 

98.8 

.29 

111 

111 


2. 12 

1. 000 

97.6 


26.8 

. 00424 

98 

18 

33 

.25 

.56 

110 

78. 0 

.38 

129 

142 


2. 19 

.908 

128 


46.0 

. 00246 

50 

22 

32 

.30 

. 69 

180 

82. 1 

.40 

176 

183 


3.37 

.962 

170 


53. 0 

.00214 

94 

19 

31 

.26 

.61 

140 

74.0 

.52 

179 

195 


3. 46 

.918 

180 


98. 2 

.00116 

92 

15 

29 

.20 

.50 

120 

62.2 

.64 

222 

262 


4. 30 

.848 

246 


50.9 

.00217 

0 

19 

28 

.29 

.67 

30 

69.6 

.22 

173 

174 

176 

3. 69 

.995 

166 


12.1 

.00914 

97 


32 




70.6 

.24 

83.9 

91 

90 

1.80 

. 923 

81.3 

No record. 

41.9 

. 00263 

.54 

is' 

29 

.27 

.61 

0 

68.3 

.36 

155 

160 

160 

3.31 

.968 

132 


51.3 

.00215 

92 

17 

27 

.26 

.62 

30 

63.5 

.47 

165 

172 

171 

3.59 

.960 

173 

Record shown in figure 4. 

116.0 

. 000953 

86 

16 

25 

.25 

.65 

0 

57.5 

.66 

235 

248 


5. 10 

.948 

260 


44. 1 

.00297 

94 

17 

33 

. 22 

.51 

0 

172 

. 67 

234 

186 

— 

4. 29 

1.258 

176 



CONSTANT LENGTH-CHORD RATIO OF 6.5— SERIES V 


1 

p 

/ sixers \ 

Percent 
F reon 

fe 

(cps) 

/ft 

(cps) . 

/. 

/« 

fe 

fR 

(deg) 

Q 

( »> ^ 

M 

V e 

(mph) 

Vr 

(mph) 

Vr' 

(mph) 

Ve 

h, . 

V, 

~Vr 

tv 

(mph) 

Remarks 


\ cu ft - / 

V sq ft / 


OoJa 


37.2 

0. 00336 

94 

30 

43 

0.31 

0.71 

10 

143 

0. 59 

197 

220 

221 

2.88 

0. 895 

201 


81.5 

. 00153 

89 

22 

40 

.23 

.55 

0 

109 

.74 

255 

318 

319 

3.73 

.804 

297 


34.7 

.00327 

96 

29 

43 

.30 

.67 


133 

.57 

193 

216 

214 

2. 78 

.893 

196 


57.4 

.00198 

95 

24 

41 

.24 

. 57 

80 

118 

.69 

234 

273 


3. 38 

.853 

252 


108 

.00105 

93 

22 

39 

.22 

.55 


90.8 

.82 

280 

363 


4.05 

.770 

345 


14.2 

. 00779 

98 

29 

37 

.37 

. 77 

0 

118 

.35 

118 

115 


2.11 

1.025 

111 

Wing failed. 

56.0 

.00197 

93 

26 

33 

.33 

.77 

0 

104 

.64 

219 

214 


3. 95 

1.023 

218 


120 

. 000923 

90 

21 

31 

.28 

.69 

0 

85.5 

.83 

291 

308 


5.24 

.945 

320 


15.8 

. 00829 

95 

39 

39 

.47 

.99 

30 

294 

.54 

181 

127 

128 

3. 11 

1.425 

113 

IModel damaged at root. 

16.7 

. 00783 

100 

39 

38 

.46 

.97 

0 

295 

.56 

186 

134 

136 

3. 05 

1.386 

122 

/Rear half separated from base. 


EFFECT MODELS— SERIES VI 


1 

P 

/slugs \ 

Percent 

Freon 

/« 

(cps) 

/ft 

(cps) 

U_ 

fn 

fe 

fu 

<P 

(deg) 

Q 

1 XI 

Ve 

(mph) 

Vr 

(mph) 

T V 
(mph) 

Ve 

btOa 

Ve 

V D 

(mph) 

Remarks 

K 

\ cu ft ) 


JR 

\ sq ft / 


1 R 











9. 15 

0. 00781 

99 

75 

76 

0. 65 

0. 89 

50 

339 

0.60 

199 

142 


2. 66 

1.40 

253 

/'Pip perpendicular to air stream. 
IModel failed. 

9. 25 

.00764 

99 

60 

78 

.51 

.70 

0 

382 

.63 

213 

146 

.... 

2.80 

1.47 

259 

(Tip perpendicular to leading edge. 
IModel failed. 

9. 55 

.00778 

99 

— 

68 

— 

— 

— 

346 

.60 

201 

127 

— 

3.02 

1.58 

229 

/Tip parallel to air stream. 
IModel failed. 

34.6 

.00205 

0 

35 

43 

.44 

.72 

— 

225 

.41 

322 

185 

189 

6.24 

1.74 

?41 

/Tip perpendicular to air stream. 
IModel failed. 

34. 1 

.00208 

0 

27 

46 

.33 

.54 

— - 

173 

.35 

278 

189 

196 

5. 21 

1.47 

348 

/Tip perpendicular to leading edge. 
IMode! failed. 

34.5 

.00207 

0 

22 

38 

.32 

.53 

0 

203 

.39 

304 

159 

159 

6. 77 

1.91 

295 

/Tip parallel to air stream. 
IModel failed. 


EFFECT OF CENTER-0 F-GRAVITY SHIFT— SERIES VII 


p 

/slugs\ 

WTft/ 

Percent 

Freon 

/« 

(cps) 

/ft 

(cps) 

/. 

fa 

fe 

fR 

<P 

(deg) 

Q 

/lb X 
\sq ft/ 

M 

V* 

(mph) 

Vr 

(mph) 

TV 

(mph) 

Ve 

bo3a 

V, 

TV 

TV 

(mph) 

Remarks 

0. 00871 

95 

12.5 

15 

0. 54 

0. 82 


153 

0. 37 

127 

231 

231 

3. 83 

0. 548 

79.9 

Model failed. 

. 00239 

0 

16 

19 

.37 

.81 

40 

109 

.28 

208 

207 

207 

3. 40 

1.000 

192 


.00177 

0 

16 

19 

.38 

.86 

20 

105 

.32 

239 

239 

239 

3. 93 

1. 000 

224 


. 0078:) 

81 

20 

21 

.47 

.94 

40 

128 

.33 

122 

120 

120 

2. 05 

1.02 

104 


. 00105 

0 

15 

18 

.35 

.83 

30 

106 

.40 

303 

308 

308 

4.97 

.985 

291 


. 00226 

0 

18 

17 

.45 

1, 09 

100 

61.5 

.20 

159 

158 

158 

2. 78 

1.01 

157 


. 00274 

76 

15 

17 

.39 

.91 

10 

58.4 

.39 

142 

141 

141 

2.54 

1.01 

139 


. 00207 

75 

14 

16 

.37 

.89 

0 

57.2 

.44 

163 

161 

161 

2. 92 

1.01 

161 


.00214 

0 

26 

36 

.42 

.72 

0 

195 

.38 

293 

415 

422 

6. 60 

.706 

245 


.00219 

0 

22 

36 

.23 

.66 

20 

151 

.33 

255 

258 

257 

3. 76 

.990 

251 


. 00224 

0 

26 

28 

.49 

.93 

20 

87. 5 

.25 

191 

176 

177 

5. 12 

1.09 

237 


.00199 

0 

26 

26 

.39 

.65 

— 

225 

.41 

324 

503 

— 

6. 73 

. 645 

267 


. 00210 

0 

23 

37 

.23 

.64 

70 

156 

.34 

264 

265 


3. 72 

.997 

257 


. 00212 

0 

23 

27 

.45 

.85 

20 

77.2 

.23 

185 

170 


5. 15 

1.09 

231 


. 00223 

0 

18 

20 

.51 

.88 

20 

61.0 

.20 

160 

160 

— 

6. 38 

1.00 

122 

] 

.00223 

0 

18 

23 

.26 

.78 

— 

62.2 

.21 

162 

139 

— 

3. 24 

1.17 

136 

{•Section reversed. 

. 00223 

0 

17 

16 

.44 

1.04 

40 

39.5 

. 17 

129 

93.2 

— 

4.78 

1.39 

110 

J 

. 00201 

0 

24 

27 

.49 

.89 

30 

258 

.44 

345 

279 

300 

5.20 

1.24 

00 

| 

. 00209 

0 

23 

26 

.48 

.86 

20 

212 

.40 

307 

186 

189 

9.15 

1.66 

00 

{■Slotted 2He inches from trailing edge. 

. 00218 

0 

23 

20 

.84 

1.03 

30 

125 

.30 

234 

121 

123 

12. 1 

1.94 

00 

J 
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